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SUMMARY 


Recorded  simulation  control  scheme  described  in  this  report  i»  an 
inexpensive  and  accurate  way  of  providing  realistic  inputs  to  a vehicle  being 
shake  tested  in  the  laboratory.  The  scheme  involves  simulating  operation 
of  the  vehicle  on  the  specified  terrain  to  generate  axle  displacement  records, 
storing  the  records  in  memory  of  a control  system,  and  using  them  to  pro- 
vide continuous  input  signals  to  the  vehicle  shaker  during  testing. 

The  computer  simulation  developed  for  the  scheme  incorporates  a 
terrain-tire- vehicle  model,  with  the  option  to  include  any  one  of  the  four 
previously  developed  tire  models;  point  contract  model,  rigid  tread  band 
model,  fixed  footprint  model  and  adaptive  footprint  model.  Developed 
separately  from  the  tire  models,  the  vehicle  model  considers  a three  axle 
military  type  truck  which  is  of  interest  to  TARADCOM.  Comparing  the 
analytical  results  generated  by  the  terrain-tire-vehicle  model  with  those 
obtained  in  field  tests  demonstrates  the  ability  of  the  model  to  simulate  the 
vehicle  motion,  particularly  when  employing  the  adaptive  footprint  or  the 
fixed  footprint  models. 

The  control  system  developed  for  storing  the  axle  displacement 
records  and  playing  them  back  incorporates  a minicomputer,  disc  memory, 
input-output  peripherals  and  digital  to  analog  converters.  In  an  illustrative 
example  described  in  the  report,  the  control  system  is  shown  to  supply 
real-time  synchronised  shaker  input  signals  from  the  digital  displacement 
records. 

Finally  the  complete  scheme  is  verified  by  comparing  the  power 
spectral  densities  of  the  inputs  to  the  vehicle  generated  by  the  scheme  with 
those  obtained  from  field  tests  on  the  vehicle. 

This  report  includes  details  of  the  model  and  the  control  system. 
However,  Instructions  for  using  the  control  system  and  the  computer  pro- 
gram are  available  in  a separate  document,  titled  "Recorded  Simulation 
Control  Scheme  - User's  Manual". 
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1.  INTRODUCTION  AND  OBJECTIVE 


In  recent  years,  laboratory  vibration  testing  of  vehicles  has  been 
widely  used  to  help  assess  structural  reliability.  In  military  service,  where 
the  dynamic  loads  are  generally  higher,  vibration  testing  has  become  a par- 
ticularly important  tool  for  reliability  evaluation,  and  has  been  used  with 
advantage  to  assess  structural  integrity,  identify  critical  stress  points,  and 
evaluate  design  improvements  (reference  1)*.  The  basic  test  procedure  con- 
sists of  shaking  a loaded  vehicle  with  hydraulic  actuators  programmed  to 
simulate  typical  field  service  conditions,  and  monitoring  the  stresses  and 
deflections  of  the  vehicle  structure.  A sketch  of  the  shaker  facility  at  the 
U.  S.  Army  Tank- Automotive  Research  and  Development  Command 
(TARADCOM)  is  shown  in  Figure  1.  The  vehicle  is  excited  by  vibratory 
inputs  applied  to  the  wheel  spindles,  as  shown.  Because  the  wheels  and  tires 
are  removed  from  the  vehicle  prior  to  shake  testing,  the  input  excitation  sig- 
nals must  correspond,  not  to  the  terrain  profile,  but  to  the  wheel  spindle 
motions  occuring  in  field  use.  There  are  two  basic  ways  in  which  these  sig- 
nals are  obtained:  either  directly  from  field  measurements,  or  from  dynamic 
computer  simulations  of  an  analytical  terrain-tire -vehicle  model.  Where 
reliable  field  data  are  readily  available,  the  former  method  can  be  followed. 

In  many  situations,  however,  such  data  are  not  easily  obtained,  and  the  latter 
method  must  then  be  used. 

There  are  many  ways  of  implementing  the  computer  controlled  Bhake 
testing  (i.  e. , using  dynamic  simulation  of  the  terrain- tire -vehicle  model). 

One  such  way,  recorded  simulation  control  scheme,  is  described  in  the 
report.  In  this  scheme,  the  vehicle  simulation  program  is  used  to  create  axle 
displacement  records  which  are  stored  in  digital  form  on  a convenient  medium, 
such  as  a magnetic  disc**.  The  records  are  then  converted  to  analog  input 
signals  and  fed  repeatedly  to  the  vehicle  shaker. 


★References  are  listed  at  the  end  of  the  report. 

★★An  intermediate  medium,  such  as  paper  tape  may  also  be  used,  before 
storing  the  data  on  magnetic  disc. 
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Before  describing  the  recorded  simulation  control  scheme,  the  imple- 
mentation options  available  to  operate  a vehicle  shaker  are  discussed  and 
advantages  of  the  selected  scheme  are  identified. 

1. 1 Implementation  Options 

There  are  several  ways  in  which  the  shake  tests  can  be  carried  out, 
depending  on  the  method  by  which  the  shaker  input  signals  are  obtained.  In 
identifying  and  evaluating  the  various  implementation  options  for  the  shaker 
input,  the  dominant  consideration  is  that  of  field  correlation,  i.  e. , the  reli- 
ability with  which  the  input  excitation  can  simulate  field  service.  Other  con- 
siderations however  are  also  important,  and  these  include  such  factors  as 
hardware  and  software  requirements,  system  complexity,  operating  conven- 
ience, etc.  The  principal  implementation  schemes  that  have  been  or  can  be 
used  to  control  the  shaker  are  summarized  in  Table  1.  They  are  described 
in  further  detail  and  evaluated  below. 

1.1.1  Direct  Field  Data  Control 


This  is  perhaps  the  simplest  of  all  the  schemes,  and  does  not 
require  any  analytical  modelling  or  computation.  It  consists  of  recording 
unsprung  mass  acceleration  (with  hub-mounted  accelerometers)  during 
vehicle  field  trials,  and  then  passing  the  signals  through  a double  integrator 
to  obtain  the  unsprung  mass  motion  (axle  motion)  signals  that  control  the 
shaker.  The  principal  advantage  of  this  method  is  that  it  is  based  on  field 
data  with  a minimum  of  signal  processing,  so  that  good  correlation  with  field 
use  can  be  expected.  However,  there  are  two  drawbacks  to  this  method 
which  are  common  to  all  schemes  using  vehicle  field  trial  data.  These  are: 

a.  Field  trial  data  are  needed  for  every  vehicle  type,  speed,  load, 
a ml  tire  pressure  of  interest. 

b>  New  concepts  and  components,  including  the  effects  of  design 
parameter  changes  cannot  be  conveniently  investigated. 
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Table  1.  Shaker  Control  Implementation  Options 


ft 

l 

l 

f. 
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Nevertheless,  this  method  will  be  Attractive  for  use  in  situations  where  teat 
data  ia  already  available,  or  where  limited  testing  of  one  or  two  vehicle 
type  a ia  contemplated. 

1«1.2  Synthesized  Field  Data  Control 

This  method  conalata  of  synthesising  the  shaker  signals  such  that 
they  correspond  statistically  to  the  field  trial  data.  It  is  best  used  when  the 
statistical  characteriatics  (PSD,  etc. } of  the  field  data  are  already  available; 
otherwise  an  initial  step  of  statistical  processing  is  required.  The  basic 
method  of  signal  synthesis  conalata  of  feeding  white  noise  into  a filter,  so  that 
tile  output  signal,  though  random,  has  the  PSD  specified  by  the  filter*.  The 
filter  characteriatics  are  chosen  so  that  the  output  PSDs  correspond  to  those 
obtained  from  field  trials.  For  statistically  independent  aignals,  this  proce- 
dure is  relatively  straightforward.  However,  for  statistically  dependent 
signals,  additional  complications  are  introduced  since  a greater  number  of 
filters  U required.  This  is  one  of  the  main  difficulties  with  this  method  for 
multi- signal  outputs:  either  a relatively  large  number  of  filters  must  be 
used,  or  statistical  dependence  (phase  relationahips)  between  the  various 
signal  pairs  cannot  be  assured.  For  signal  pairs  whose  statistical  dependence 
is  known  to  be  weak,  the  numbor  of  filters  can  be  minimized  by  considering 
the  signals  independent.  However,  since  the  degree  of  dependence  changes 
with  the  terrain  and  vehicle  characteristics,  speed,  load,  etc,,  significant 
statistical  analysis  will  be  needed  to  identify,  for  esch  case,  the  outputs  that 
can  be  synthesized  independently. 


| 


! 


i 


This  implementation  scheme  is  attractive  for  use  primarily  in 
situations  where  a relatively  small  number  of  signals  are  to  be  synthesised, 
and  where  the  statistical  characteristics  of  field  trial  data  are  readily  avail- 
able, It  will  also  be  attractive  for  use  when  existing  hardware  is  suitable. 


♦The  filters  and  noise  sources  can  be  implemented  with  analog  hardware  or 
with  digital  software,  as  shown  in  Table  1. 
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1.1.3  Direct  Simulation  Control 


Direct  simulation  control  is  a method  by  which  the  shaker  is 
controlled  on-line  by  a computer.  It  is  a scheme  which  requires  no  field 
trial  data,  but  generates  the  input  signals  by  means  of  a dynamic  terrain- 
tire-vehicle  simulation.  Its  implementation  is  relatively  straightforward, 
but  it  requires  use  of  an  adequate  on-line  computer  during  the  entire  shake 
test.  For  facilities  with  access  to  an  on-line  computer,  this  method  can  be 
very  attractive  because  of  its  simplicity  of  implementation.  Where  an  ade- 
quate on-line  simulation  cannot  be  conveniently  performed,  a modified 
scheme  which  retains  the  basic  simplicity  can  be  used.  This  scheme  is 
described  in  the  following  section. 

1.1.4  Recorded  Simulation  Control 

This  scheme  carries  out  the  dynamic  vehicle  simulation  similar 
to  the  direct  simulation  control  method,  but  records  the  output  on  a multi- 
track data  recorder  (e.g. , on  magnetic  disc  or  tape).  The  records  are  then 
played  back  to  control  the  shaker.  Thus  by  using  the  recorder,  the  need  for 
an  on-line  computer  is  eliminated.  This  improvement  makes  the  recorded 
simulation  control  scheme  particularly  attractive  for  use.  The  hardware 
requirements  are  modest,  since  the  simulation  can  now  be  carried  out  at  a 
batch  processing  facility.  Field  trial  data  io  not  needed,  and  a variety  of 
vehicles,  types  and  loading  conditions  can  be  conveniently  tested  using  the 
simulation.  No  statistical  processing  is  required,  and  since  the  shaker  sig- 
nals are  generated  directly  by  a time  simulation,  the  statistical  dependence 
of  all  the  signals  will  be  preserved.  For  the  above  reasons,  the  recorded 
simulation  control  scheme  is  well  suited  to  situations  where  different  types 
of  vehicles  (including,  specifically,  multi-axle  vehicles)  have  to  be  tested 
under  various  load  and  speed  conditions,  and  where  parameter  and  component 
changes,  including  new  concepts,  are  to  be  evaluated. 
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1.1.5  Synthesised  Simulation  Control 


The  synthesised  simulation  control  scheme  is  essentially  simi- 
lar to  the  synthesised  field  data  control  scheme  discussed  earlier;  the  only 
difference  being  that  the  statistical  characteristics  of  wheel  motion  are 
obtained  by  analytical  simulation  rather  than  from  field  trial  data.  The  basic 
method  of  synthesising  the  signal  by  feeding  white  noise  through  a filter 
remains  unchanged.  As  before,  hardware  requirements  (filters  and  noise 
generators,  or  a hybrid  computer)  can  be  high  if  statistical  dependence  is  to 
l be  preserved.  Although  field  trial  data  will  not  be  needed,  significant  com- 

[ putation  and/or  statistical  processing  may  be  required,  particularly  for 

multi-axle  vehicles.  However,  the  convenience  of  "hardwired"  electronics 

£ 

or  computer  software*  make  it  particularly  attractive  for  long-term  testing 
of  a limited  number  of  vehicle  types  and  load  conditions. 

L 

!' 

1.1.6  System  Evaluation 

i 

Evaluation  and  selection  of  recommended  schemes  for  implemen- 
tation of  TARADCOM  must  be  based  on  a number  of  factors,  including: 

The  ability  of  the  scheme  to  test  vehicles  for  which  field  trial  data 
is  unavailable. 

The  capability  of  evaluating  a variety  of  vehicle  types,  including 
specifically  multi-axle  trucks  with  a relatively  large  number  of 
wheel  excitation  inputs. 

The  hardware  and  software  requirements,  and  their  availability 
at  TARADCOM. 

System  complexity,  ease  of  implementation,  operating  convenience, 
etc. 


♦The  noise  generator  and  filter  system  where  modeling  changes  are  made  by 
adjusting  dial  settings  or  by  changing  coefficients  in  a digital  computer 
program. 
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The  first  two  scheme • (direct  field  data  control  and  synthesised 
field  data  control)  described  in  the  previous  section  require  vehicle  field 
trial  data,  and  their  use  at  TARADCOM  depends  upon  the  availability  of  such 
data.  The  third  scheme  (direct  simulation  control),  while  straightforward 
in  implementation,  requires  an  adequate  on-line  analog  (or  hybrid)  com- 
puter • the  availability  of  which  will  be  a primary  factor  in  deciding  upon 
the  use  of  this  scheme  at  TARADCOM.  The  fifth  scheme,  the  synthesised 
simulation  control  scheme,  has  the  advantage  of  not  requiring  test  data,  and 
will  be  convenient  to  use  for  long-term  testing  of  a limited  number  of  vehicle 
types.  However,  for  multi-axle  trucks,  implementation  will  be  relatively 
complex  and  hardware  requirements  may  be  high. 

The  recorded  simulation  control  scheme  has  many  desirable 
characteristics  which  make  it  well  suited  for  implementation.  It  does  not 
require  field  trial  data.  Computation  requirements  are  minimised,  since 
no  statistical  data  processing  is  necessary.  The  dynamic  vehicle  simulation 
can  be  carried  out  off-line  at  a batch-processing  facility,  eliminating  the 
need  for  an  on-line  computer.  The  system  can  conveniently  generate  multiple 
signal  outputs  and,  unlike  signal  synthesis  schemts,  no  special  precautions 
are  necessary  to  preserve  tho  phase  relationships  of  the  various  wheel  exci- 
tation inputs.  The  hardware  requirements  needed  to  implement  the  scheme 
are  modest  - the  only  major  item  required  being  a multi-track  data  storage 
and  retrieval  system.  Therefore,  this  scheme  is  selected  for  implementa- 
tion and  the  rest  of  the  report  deals  with  the  details  of  the  scheme  in  terms 
of  the  dynamic  vehicle  simulation  program  and  a digital  controller  which 
performs  the  data  storage  and  retrieval  functions. 


1 . 2 The  Recorded  Simulation  Control  Scheme 

There  are  several  ways  of  implementing  the  recorded  simulation  con- 
trol scheme,  with  the  same  basic  approach  of  recording  the  computer 
generated  axle  displacement  signals  and  playing  them  back  at  the  time  of 
shake  test.  The  implementation  scheme  selected  in  this  effort  is  shown  in 
Figure  2. 


[ ! 

'■< 

As  shown  in  the  figure,  the  axle  motion  records  are  generated  by  a 
dynamic  vehicle  simulation  computer  program.  Data  necessary  to  execute 

I the  program  includes 

( { 

I a.  Vehicle  data  (mass,  GG  location,  wheel  base,  etc. ) 

| b.  Operation  data  (terrain  profile  and  speed  of  operation), 

: 

, » 

! 

The  axle  displacement  records  generated  by  the  simulation  are  stored 
on  paper  tapes  and  transferred  to  magnetic  disc  of  a digital  controller 
through  adequate  software.  The  data  stored  on  the  magnetic  disc  are  re- 
trieved  at  the  time  of  running  the  shake  test  and  the  digital  signals  are 
converted  to  analog  signals  using  a bank  of  digital- to-analog  signal  converters. 
The  analog  signals  are  then  fed  to  the  shaker  through  appropriate  filters  and 
amplifiers. 

The  record  stored  on  the  disc  typically  lasts  less  that  a minute  in 
real-time:  therefore,  the  process  is  repeated  until  the  testing  is  completed. 
Also,  the  simulation  program  is  executed  each  time  the  vehicle  test 
parameters  (terrain  profile,  speed,  load)  are  changed,  and  the  axle  displace- 
ment records  on  the  disc  are  updated  to  show  the  changes. 


*j 


1.3  Objective  and  Organisation  of  the  Report 

Objective  of  the  report  is  to  introduce  the  recorded  simulation  control 
scheme,  describe  it  and  demonstrate  its  effectiveness  through  an  illustrative 
example.  The  report  describes  the  scheme  in  terms  of  its  two  steps*  (a)  gen- 
erating axle  displacement  records  using  the  vehicle  simulation  program,  and 
(b)  obtaining  input  signals  for  the  vehicle  shaker  from  the  axle  displacement 
records.  The  vehicle  simulation  program  incorporates  a terrain-tire-vehicle 
model  developed  in  Chapter  Z,  There  are  three  components  in  this  modelt 
terrain,  tire  and  vehicle.  The  terrain  model  generates  a terrain  profile  from 
information  on  its  spectral  density.  The  vehicle  model  developed  in  this 
report  allows  simulation  of  vehicles  similar  to  the  5-ton  M-809  cargo  truck. 
Tire  represents  an  interface  between  the  two  models  and  it  is  modelled  in 
four  different  ways  (see  reference  2)j 
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a.  Point  contact  model 

b.  Rigid  tread  band  model 
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c.  Fixed  footprint  model 

d.  Adaptive  footprint  model. 

Each  component  of  the  overall  ter  rain- tire -vehicle  model  is  developed 
separately  ao  that  it  can  be  used  independently  in  different  applications. 


Chapter  2 also  discussed  results  obtained  from  a typical  simulation 
run  using  a computer  programs  incorporating  the  model  developed.  The 
simulation  results  are  then  compared  with  corresponding  experimental 
results  to  formulate  guidelines  for  selecting  a tire  model  and  verify  the 
overall  model. 

The  controller  needed  to  store  and  retrieve  the  axle  displacement 
records  generated  by  tho  vehicle  simulation  is  described  in  Chapter  3.  The 
necessary  hardware  components  and  their  functions  are  discussed  along  with 
key  steps  in  operation  of  the  controller.  Details  of  the  controller  operation 
are.  however,  not  discussed  because  they  are  available  in  the  User's  Manual 
(reference  4). 

Chapter  4 deals  with  verification  of  the  complete  recorded  simulation 
control  scheme  through  comparison  of  the  spectral  densities  of  the  axle 
displacement  signal  with  those  obtained  from  field  data.  Finally,  conclusions 
and  recommendations  are  discussed  in  Chapters  £ and  6. 


♦The  computer  program  itself  is  described  in  a separate  document 
(reference  4). 
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2.  SIMULATION 


Development  of  a reliable  model  for  simulating  the  vehicle  performance 
was  initiated  in  the  first  phase  of  the  project  (reference  2),  It  was  realiaed 
early  in  the  effort  that  a good  model  of  tire  is  necessary  for  accurate  per- 
formance prediction  and  therefore  in  the  first  phase  four  different  tire  models 
were  developed.  These  models  were  incorporated  in  a simple  two  degree  of 
freedom  model  and  initial  simulation  results  were  obtained  to  provide  insight 
into  tire  model  behavior. 

These  tire  models  are  now  incorporated  in  a six  degree  of  freedom 
vehicle  model  which  can  simulate  the  heave-pitch  motion  of  a three  axle  mili- 
tary type  truck.  The  vehiole  has  been  selected  so  that  the  results  obtained 
from  the  simulation  can  be  compared  with  the  available  experimental  results. 
Also,  the  input  signals  generated  by  the  scheme  can  then  directly  be  used  to 
test  similar  vehiclss  on  the  six  actuator  shaker  available  at  TARA  DCOM. 
Because  of  the  non-gene raliaed  nature  of  the  vehicle  model,  modifications 
are  necessary  for  testing  other  types  of  vehicles  using  the  scheme  described 
in  the  report.  However,  tin  terrain  and  the  tire  models  are  developed  so 
that  they  can  be  used  independently  in  any  application. 

The  simulation  results  described  in  the  report  are  obtained  using  a 
computer  program  which  incorporates  the  model.  Details  of  the  computer 
program  are  available  in  the  User’s  Manual  (reference  4)  and  therefore  are 
not  repeated  here. 

2. 1 Model  Development 

The  tire,  vehicle  and  terrain  models  are  summarised  in  the  sections 
that  follow.  Additional  details,  including  equations,  for  the  tire  and  vehicle 
models  are  presented  in  the  Appendices, 
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2.1.1  The  Tire  Model 


Four  bailie  tire  model*  were  formulated  in  Phaae  1.  ranging  in 
sophistication  from  the  aim  pie  point  contact  model  to  the  more  advanced 
adaptive  footprint  model.  These  model*  were  developed  no  that  the  moat 
economical  tire  simulation  can  be  selected  for  the  task  at  hand.  The  models 
are  summarised  in  Table  2 and  described  briefly  below.  Further  details 
can  be  found  In  reference  2, 

2.  1.  1. 1 Thu  Point  Contact  Tire  Model 


The  point  contact  model  is  represented  by  a parallel 
spring-dashpot  combination  that  transmits  the  support  force  from  the  terrsin 
to  the  vehicle  and  contacts  the  ground  through  a point  follower.  The  spring 
stiffness  is  chosen  to  simulate  the  effects  of  inflation  pressure  and  carcass 
elasticity.  The  damping  provides  the  energy  dissipation  caused  by  carcass 
deformations.  Terrain  contact  occurs  at  a single  point  vertically  beneath 
the  wheel  center.  Dynamic  support  forces  occur  due  to  deflection  of  the 
spring  and  daahpot  caused  by  motion  of  the  whoel  spindle  relative  to  the 
terrain.  Fore-and-aft  forces  are  obtained  from  this  model  by  assuming  that 
the  resultant  tire  force  is  always  normal  to  the  local  terrain  surface.  When- 
ever the  profile  surface  is  inclined  to  the  horiaontal,  a fore-and-aft  force  is 
generated  which  is  related  to  the  vertical  force  through  the  tangent  of  the 
local  surface  angle.  The  tire  mass  is  concentrated  at  the  wheel  center,  and 
the  terrain  follower  is  free  to  leave  the  ground  to  simulate  wheel  hop. 

2, 1.1.2  Riaid  Tread  Band  Tire  Model 

The  rigid  tread  band  tire  model  is  similar  to  the 
point  contact  model,  except  that  the  terrain  contact  point  is  not  constrained 
to  lie  vertically  beneath  the  wheel  center  but  is  free  to  move  forward  or  aft 
depending  on  the  local  profile  slope.  This  model  is  equivalent  to  a point 
contact  tire  with  the  terrain  contact  point  located  at  the  center  of  a rigid 
rolling  wheel  or  circular  tread  band.  In  this  case,  the  wheel  center  sees  a 
motion  that  Is  in  general  different  from  the  terrsin  profile,  due  to  filtering 
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effects  of  the  rolling  wheel.  On  rough  ground,  the  filtering  causes  a high 
Attenuation  of  profile  irregularities.  For  gradual  changes  in  profile  eleva- 
tion, filtering  becomes  insignificant,  and  the  two  models  give  equivalent 
results. 


As  with  the  point  contact  model,  tha  support  force  is 
transmitted  to  the  vehicle  through  the  parallel  spring -dashpot  combination 
that  simulates  the  inflation  pressure  and  carcass  forces.  Fore-and-aft 
forces  are  obtained  by  assuming  that  the  total  tire  force  always  acts  normal 
to  the  terrain  at  the  contact  point.  The  tire  mass  is  concentrated  at  the 
wheel  center,  and  the  tread  band  is  free  to  leave  the  ground  to  simulate 
wheel- hop. 


Z. 1,1. 3 Fixed  Footprint  Tire  Model 

The  fixed  footprint  tire  model  interacts  with  the  ter- 
rain through  a footprint  of  constant  sise,  independent  of  the  tire  deflection. 
Inflation  pressure  and  carcass  forces  are  simulated  by  the  spring  and  damper 
elements  distributed  uniformly  over  the  contact  length.  The  finite  footprint 
area  provides  this  model  with  the  ability  to  envelop  small  irregularities 
through  local  deformations  within  the  footprint.  In  fact,  * this  model  is 
equivalent  to  a point  contact  model  in  which  the  local  terrain  elevation  is 
replaced  by  the  average  terrain  elevation  taken  across  the  footprint  length. 
Like  the  rigid  tread  band  model,  the  fixed  footprint  model  also  filters  the 
terrain  profile,  except  that  the  filter  characteristics  are  governed  by  the 
footprint  length  rather  than  the  tire  radius.  Unlike  the  rigid  tread  band  model, 
however,  the  footprint  center  is  constrained  to  lie  vertically  below  the  wheel 
centert  independent  of  the  local  profile  slope. 

The  fore-and-aft  tire  force  is  obtained  by  assuming 
that  the  footprint  force  acts  perpendicular  to  the  local  profile  and  the  foot- 
print is  free  to  leave  the  ground  to  simulate  wheel -hop. 


♦With  linear  stiffness  and  damping. 
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2. 1.1.4  Adaptive  Footprint  Tire  Model 


The  adaptive  footprint  tire  model  consist*  of  a flexible 
tread  band  inflated  by  internal  pressure  and  linked  to  the  wheel  center  by 
angularly  distributed  stiffness  and  damping  which  simulate  carcass  and  tread 
stiffness.  As  the  wheel  passes  over  terrain  irregularities,  the  tread  in  the 
footprint  region  deforms,  and  gives  rise  to  the  pneumatic  and  carcass  com- 
ponents of  tire  force.  The  former  is  given  by  the  inflation  pressure  acting 
over  the  contact  area.  The  latter  is  found  from  the  resultant  of  the  spring 
and  damping  forces  in  the  footprint.  In  general,  the  resultant  tire  force  will 
not  be  vertical  because  of  the  existence  of  a non- planar  footprint,  and  the 
calculation  of  the  appropriate  components  will  allow  determination  of  the 
vertical  and  fore-and-aft  tire  force. 

Like  the  fixed  footprint  model,  this  model  has  the 
ability  to  envelop  small  irregularities  through  local  deformations  within  the 
footprint.  However,  the  key  feature  of  this  model  is  that  the  footprint  size 
and  orientation  relative  to  the  wheel  center  changes  depending  on  the  tire 
deflection  and  terrain  profile.  As  with  the  other  models,  the  tire  mass  is 
concentrated  at  the  wheel  center  and  the  tire  is  free  to  leave  the  ground  to 
simulate  wheel  hop. 

2.1.2  The  Vehicle  Model 


The  vehicle  model  has  been  formulated  to  allow  simulation  of  a 
variety  of  trucks  of  interest  to  the  military.  The  baseline  vehicle  which 
guided  simulation  development  and  was  used  to  validate  the  model  was  the 
M-809  5 ton,  three  axle,  6x6  cargo  truck.  This  vehicle  was  selected  both 
because  field  trial  data  were  available,  and  because  the  basic  vehicle  con- 
figuration is  representative  of  a variety  of  multi-axle  military  trucks. 

A sketch  of  the  vehicle  is  shown  in  Figure  3a.  The  front  wheels 
are  attached  to  the  frame  through  a semi-eliptical  leaf  spring  and  dashpot 
suspension.  The  rear  suspension  (Figure  3b)  consists  of  a pivoted  leaf  spring 
assembly  (bogie)  attached  to  the  chassis  through  a central  shaft  (pivot)  and 
connected  at  t^e  spring  ends  to  the  middle  and  rear  wheel  spindles.  Bounce 


16 


motion  of  the  rear  wheels  thus  results  in  two  types  of  suspension  motion: 
flexure  of  the  leaf  spring,  and  rotation  of  the  bogie  about  the  pivot.  For 
symmetrical  vibration  - i.  e. , when  the  rear  wheel  pair  vibrate  in  phase  - 
the  spring  is  flexed  without  bogie  rotation.  For  antieymmetrical  vibration  - 
i.  e. , when  the  rear  wheel  pair  vibrate  out  of  phase  - the  bogie  rotates  without 
spring  flexure.  In  the  real  case,  when  the  vehicle  is  going  over  a bumpy  road, 
both  symmetrical  and  antisymmetrical  components  will  exist,  and  spring 
flexure  and  bogie  rotation  will  occur  simultaneously.  However,  excessive 
jounce  motion  of  the  rear  wheels  is  prevented  by  stops  attached  to  the  main 
frame. 


The  analytical  model  for  the  truck  is  shown  in  Figure  3c.  The 
vehicle  body  is  represented  by  a rigid  beam  of  mass  M^  and  rotational  inertia 
1^.  The  unsprung  mass  of  the  front,  middle  and  rear  wheels  (M^,  M3,  M^) 
consists  of  the  effective  mass  of  the  tire(s)  and  the  mass  of  the  wheel(a), 
sunpeusion  and  other  components  that  move  with  the  wheel.  The  vertical  tire 
force  acts  directly  on  the  unsprung  mass.  The  fore-and-sft  tire  force  is 
transmitted  to  the  vehicle  body  through  suspension  elements  not  shown  in  the 
figure.  The  front  suspension  leaf  spring  is  modelled  as  a linear  spring  (K2) 
with  a hardening  characteristic  for  large  motions  (elastic  stops),  as  shown 
in  Figure  4a,  When  the  spring  excursion  exceeds  the  jounce  or  rebound 
clearance,  its  stiffness  increases  by  a factor  of  n (typically,  n = 10)  due  to 
stop  contact.  Dry  friction  damping  (B2j)  is  also  included  in  the  leaf  spring 
model,  as  shown  in  Figure  4b.  The  dry  friction  force  has  a ralue  of  5 per- 
cent of  the  instantaneous  spring  force  (before  stop  contact),  and  opposes  the 
velocity.  The  front  shock  absorber  is  modelled  as  a daahpot  with  damping 
coefficient  B2y,  in  jounce  and  B2yr  in  rebound,  as  shown  in  Figure  4c.  The 
rear  suspension  bogie  assembly  is  modelled  as  a rigid  beam  (mass  M^, 
rotational  inertia  1^)  pivoted  at  the  center,  with  two  springs  (K^)  attached 
to  the  extremities.  The  beam  1 'presents  the  kinematic  configuration  and 
rotational  inertia  of  the  bogie  during  antisymmetrical  wheel  vibration  (i,  e. , 
for  bogie  rotation).  The  spring  elements  (K^)  represent  each  half  of  the 
rear  leaf  spring.  As  with  the  front  spring,  the  rear  spring  is  also  modelled 
as  a linear  spring  with  elastic  stops  and  dry  friction  damping  (Figures  4a 
and  4b).  No  viscous  damping  is  included  in  the  rear  suspension  model, 
because  of  the  absence  of  a shock  absorber.  Frame  mounted  stops  which 
prevent  excessive  rear  wheel  jounce,  and  tire  stops  that  limit  the  maximum 
tire  deflection  are  also  included  in  the  model.  The  force-deflection 
characteristics  of  these  elements  are  shown  in  Figure  S. 


18 


Force 


5 percent  of 

1 Instantaneous  Spring 

Force 

Velocity 

s 

v—  5 percent  of 
Instantaneous 
Spring  Force 

Force 


Slope  - Damping 

Constant,  B-  , 
2vj 


. Slope  - Damping 
Constant,  B2yr 


(b)  Dry  Friction  Damping 


(c)  Shock  Absorber  Damping 


Figure  4,  Force  Characteristics  of  Suspension 
Spring  and  Damping  Models 
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Figure  5.  Force  Characterietic*  of  Frame  Stop 
and  Tire  Model* 


The  vehicle  model  described  above  has  been  developed  primar- 
ily for  the  initial  bounce-pitch  simulation!  needed  to  verify  the  tire  models 
and  carry  out  the  demonstration  shake  testa.  Subsequently,  the  roll  mode 
can  be  included  by  providing  separate  left  and  right  wheel  inputs  and  including 
vehicle  body  roll  inertia. 

2.1.3  The  Terrain  Model 

The  terrain  model  developed  for  the  simulation  consists  of  a 
series  of  profile  elevation  coordinates  whose  values  are  determined  statis- 
tically, and  can  represent  several  types  of  unimproved  roads  of  varying 
degrees  of  roughness  (reference  2),  This  terrain  model  has  been  chosen 
primarily  to  make  the  verification  runs  of  the  vehicle  simulation  program, 
and  for  the  initial  shake  tests.  Subsequently,  the  analytically  obtained  ter- 
rain coordinates  can  be  replaced  conveniently  with  real  terrain  data. 

A 500  foot  length  of  profile  has  been  generated  and  used  for  the 
model  validation  runs.  An  rms  roughness  of  1 inch  was  chosen  for  the  simu- 
lated profile,  so  that  its  power  spectral  density  (PSD)  generally  agrees  with 
Aberdeen  Proving  Grounds  test  track  data  (reference  3),  for  which  the  M-809 
truck  trial  results  were  available.  The  statistical  terrain  characteristics 
are  summarised  in  Table  3,  The  PSD  of  the  profile  is  shown  in  Figure  6. 

As  shown  in  the  figure,  the  PSD  is  plotted  as  a function  of  a 
reduced  frequency  ft  (=  ?tt/A),  where  X is  the  wavelength  of  the  irregularity 
component.  This  frequency  can  be  converted  to  frequency  in  Hx  by  multiplying 
by  V/2ti  . Corresponding  PSD  is  obtained  by  multiplying  the  PSD  shown  by  2ir/V. 

The  straight  line  on  the  figure  represents  the  analytically  specified 
PSD.  The  points  were  found  by  computing  the  PSD  from  the  actual  elevation 
coordinates.  The  low  frequency  scatter  is  due  in  part  to  the  reduced  number  of 
complete  cycles  available  for  calculating  the  mean  square  values,  and  can  be 
reduced  by  Increasing  the  profile  length.  The  dynamic  simulation  results  and 
experimental  comparison  are  presented  in  the  sections  that  follows. 


Table  3.  Statistical  Terrain  Characteristics 


rms  Roughness,  yrmB 

1 inch 

PSD  Relationship,  SyQ  (SI) 

A/S32 

Amplitude  Probability  Distribution 

Gaussian 

Upper  Cutoff  Wavelength,  Xj 

57  feet 

Lower  Cutoff  Wavelength,  X2 

0. 177  feet 

Vehicle  Speed,  V 

18  mph 

Excitation  Frequency  Range 

0.46  - 150  Ha 

2.2  Simulation  Results 


The  terrain,  tire  and  vehicle  models  described  previously  were  inte- 
grated into  a dynamic  bounce-pitch  simulation  suitable  for  obtaining  the 
shaker  input  signals.  Three  programs  have  been  developed:  program  GRND 
which  generates  the  terrain  profile,  program  TIRE  which  determines  the 
vehicle  motion  time  history,  and  program  SPEC  which  computes  the  PSDs  of 
vehicle  motion.  These  programs,  and  their  method  of  use,  is  described  in 
detail  in  a separate  document  (reference  4).  A test  case  simulation  using 
these  programs  was  carried  out  to  verify  the  models  and  provide  the  neces- 
sary data  for  the  initial  demonstration  shake  tests.  The  simulations  were 
carried  out  for  an  M-809  truck  moving  at  18  mph  over  a 500  foot  long  track 
of  1 inch  rms  roughness.  The  parameters  of  the  tire  and  the  vehicle  used  for 
the  simulation  are  summarised  in  Table  4,  Table  5 shows  additional  tire 
parameters  which  are  required  for  the  more  sophisticated  tire  models. 

A short  extract  of  the  time  history  of  tire  force  obtained  with  each  of 
the  tire  models  is  shown  in  Figure  7.  Examination  of  the  force  response  of 
the  point  contact  tire  shows  that  it  is  significantly  higher  than  that  of  the 
other  models,  and  has  a greater  high  frequency  content.  Another  character- 
istic seen  from  the  time  simulation  is  that  wheel  hop  (aero  tire  force)  is 
important.  In  the  10  foot  section  of  record  shown,  the  point  contact  model 
spends  almost  60  percent  of  the  time  out  of  ground  contact.  The  other  two 
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models  (rigid  tread  band  and  fixed  footprint)  also  show  wheel  hop,  although 
in  decreasing  amounts.  Although  the  adaptive  footprint  simulation  is  free  of 
wheel  hop  for  the  10  foot  section  shown,  the  wheel  hop  condition  is  closely 
approached  on  two  occasions.  The  reason  for  the  reduction  in  wheel  hop  with 
the  more  advanced  models  is  that  they  filter  the  short  wavelength  irregular!* 
ties,  and  hence  attenuate  the  high  frequency  excitation  which  causes  wheel 
hop.  The  filtered  ground  profiles  for  the  rigid  tread  band  and  fixed  footprint 
tire  models  are  shown  in  Figure  8 along  with  the  actual  profile  for  compari- 
son. The  reduction  in  high  frequency  content,  especially  for  the  fixed  foot- 
print model  is  readily  apparent. 


Table  4.  Summary  of  Vehicle  Parameters 


Vehicle  Type  - M-809,  5 ton,  6x6  cargo  truck 

Tire  - Goodyear  11.00-18,  Military  NDCC 

Vehicle  Body  Parameters 

Symbol 

Value 

Unit 

Mass  of  Hull 

Mh 

718.6 

Slugs 

Mass  of  Bogie 

“b 

6.0 

Slugs 

Front  Unsprung  Mass  (each  side) 

m2 

40.6 

Slugs 

Middle  Unsprung  Mass  (each  side) 

m3 

42.1 

Slugs 

Bear  Unsprung  Mass  (each  side) 

m4 

43.6 

Slugs 

Pitch  Inertia  of  Hull  (about  eg) 

34924 

Slug  *ft2 

Pitch  Inertia  of  Bogie  (about  pivot) 

*b 

3.0 

Slug,  ft2 

Front  Suspension  to  Hull  CG  Length 

L2 

12.  5 

ft 

Bogie  Pivot  to  Lull  CG  Length 

L3 

2.4 

ft 

Hull  CG  Height 

GG 

0 

ft 

Offset,  Bogie  Pivot  to  Hull  Plane 

S34 

0 

ft 

Offset,  Bogie  Pivot  to  Bogie  CG 

Sb 

0 

ft 

Bogie  Pivot  to  Front  Stop  Length 

sf. 

2.2 

ft 

Bogie  Pivot  to  Bear  Stop  Length 

Sb. 

2.2 

ft 

Half  Bogie  Length 

a 

2.2 

ft 
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Table  4.  Summary  of  Vehicle  Parameter*  (Coat) 


Vehicle  Body  Parameter*  (Cont) 

Symbol 

Value 

Unit 

Stop  Length 

0.25 

ft 

Stop  Clearance  (Unloaded) 

St  ^ 
JL*P 

0.39 

ft 

Suspension  Parameters 

Symbol 

Value 

Unit 

Front  Suspension  (each) 

Stiffness 

K* 

29976 

lb/ft 

Viscous  Damping  Coefficient  (Jounce) 

B2vj 

185 

lb*  eec/ft 

Viscous  Damping  Coefficient  (Rebound) 

B2vr 

600 

lb'  aac/ft 

Jounce  Travel  (Unloaded) 

Jn2 

0.5 

ft 

Rebound  Travel  (Unloaded) 

Rb2 

0.5 

ft 

Dry  Friction  Coefficient 

B2f 

e 

Middle  and  Rear  Suspension  (each) 

Stiffness 

K3 

87000 

lb/ ft 

Jounce  Travel  (Unloaded) 

Jn3 

0.7 

ft 

Rebound  Travel  (Unloaded) 

Rb3 

0.7 

ft 

Dry  Friction  Coefficient 

S3f 

* 

Tire  Parameters 

Symbol 

Value 

Unit 

Stiffness 

k 

39000 

lb/ ft 

Damping 

b 

31 

lb*  sec/ft 

Radius 

r 

1.67 

ft 

Contact  Length 

L 

1.03 

ft 

Effective  Width 

B 

0.52 

ft 

Inflation  Pressure 

Pi 

30 

psi 

Tire  Deflection  Limit 

0.4 

ft 

^Selected  each  that  the  dry  friction  damping  force  i*  equal  to  5 percent  of 
the  spring  force,  with  direction  opposing  the  velocity. 
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Table  5.  Model  Parameter*  for  11.00-18  Tire 
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Point  Contact  Modal 


Rigid  Tread 
Band  Model 


X Diapln 

X Diapln 

Fixed  Footprint  Model 

2.4 

Adaptive  Footprint  Model 

454  456  458  460 

X Diapln  X Diapln 

Notes  X axie  - diatance,  ft 

Y axis  - force/equilibrium  force 

(a)  Vertical  Forco 

Figure  7.  Time  Hiatory  of  Tire  Force 
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Point  Contact  Modal 


Rigid  Tread  Band  Model 


X Diapln 

X Diapln 

0.  28 

Fixed  Footprint  Model 

_ Adaptive  Footprint  Model 

X Diapln  X DUpln 

(b)  Fore-and-aft  Force  Note:  X axia  - diatance,  ft 

Y axis  - force /equilibrium  force 

Figure  7.  Time  Hiatory  of'  Tire  Force  (Cont) 
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3 


X Diepln 


Note: 

X axis 
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distance,  ft. 
elevation,  ft, 


Figure  8,  Actual  and  Filtered  Ground  Profile 


Additional  conclusions  regarding  the  tire  model*  can  be  drawn 
from  the  FSDi  of  the  time  functions.  The  PSDs  of  tire  force  are  shown 
in  Figures  9 and  10.  Figure  9 shows  the  PSD  of  vertical  tire  force 
obtained  with  each  tire  model.  There  are  three  distinct  peaks  that  can 
be  observed.  The  first  peak,  at  about  0.5  rad/ft  corresponds  to  the 
vehicle  body  bounce/pitch  natural  frequency.*  The  second  peak,  occur- 
ing  at  about  1.5  rad/ft,  corresponds  to  the  bogie  pitch  natural  frequency 
caused  by  antisymmetric  (out -of -phase)  vibration  of  the  two  rear  axles. 
The  third  peak,  occuring  at  about  2.5  rad/ft,  corresponds  to  the  conven- 
tional unsprung  mass  natural  frequency. 

As  an  initial  check  of  the  program  formulation,  the  above  natural 
mode  frequencies  have  been  compared  with  those  found  from  a simplified 
approximate  analysis.  The  approximate  expressions  for  natural  frequency 
(u^)  are  obtained  by  choosing  the  appropriate  stiffness  (K)  and  inertia  (M) 
and  then  determining  un  a J K/M  as  follows. 


(1) 


,<2> 


(3) 


Body  Bounce, 


Sab 


J l Kzk  4ISk  \ / 

*\z\ K-i*  + sjrar)  / “h 


Body  Pitch, 


Bogie  Pitch, 


“np  = V2WM4 


Since  body  bounce  and  pitch  natural  frequencies  for  the  M809  truck  are 
very  close  together,  two  separate  peaks  are  not  always  observed. 
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(4) 


Front  Wheel  Bounce,  u^bf  — V(^2+k)/M-> 

Rear  Wheel  Bounce,  w'nbr  — V (K^+ZkJ/M^  (5) 

. 

These  expressions,  evaluated  from  Table  4,  are  summarized  in 

Table  6.  As  can  be  seen,  the  simulation  results  are  in  close  agreement 

* * 

with  the  approximate  analysis,  which  serves  as  an  initial' check  of  the 
program  formulation. 

The  simulation  results  (Figures  9 and  10)  show  that  the  simpler 
tire  models  generally  overestimate  the  transmitted  tire  force.  Although 
at  low  frequencies  all  models  give  similar  results,  the  differences 
become  apparent  in  the  middle  and  upper  frequency  ranges  of  interest. 

For  instance,  point  contact  model  predictions  at  the  higher  frequencies 
can  be  up  to  an  order  of  magnitude  greater  than  those  of  the  more 
advanced  models.  The  rigid  tread  band  model  also  gives  generally 
higher  force  estimates,  specially  for  the  fore-and-aft  force.  The  fixed 
footprint  and  adaptive  footprint  simulations,  however,  are  in  good  agree- 
ment with  each  other,  and  give  the  lowest  force  predictions,  because  of 
their  ability  to  envelop  small  irregularities  through  local  footprint 
deformations. 

Similar  conclusions  regarding  the  various  tire  models  can  be 
reached  by  examining  the  other  simulation  results  (CG  acceleration, 
wheel  acceleration,  etc.).  These  results  confirm  that  the  adaptive  foot- 
print and  fixed  footprint  tire  models  agree  closely,  and  predict  lower 
accelerations  than  the  rigid  tread  band  and  point  contact  models.  For 
reference,  all  the  principal  simulation  results  are  presented  in  Appendix  C. 

2.3  Simulation  Verification 


The  basic  objective  in  developing  the  present  terrain-tire-vehicle 
simulation  is  to  provide  a means  of  generating  realistic  shake  test  input 
signals.  Therefore,  before  carrying  out  the  tests,  it  is  necessary  to 
verify  the  model  oredictions  with  field  trial  data  so  that  the  shake  test 
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Table  6.  Comparison  of  Simulation  Results  vith 
Approximate  Analysis 


Natural  Mode 

Natural  Frequency  (He) 

Bounce -Pitch 
Computer 
Simulation 

Approximate 

Analysis 

Body  Pitch, 

np 

2.2 

2. 1 

Body  Bounce, 

2.2 

2.6 

Bogie  Pitch,  w'np 

6.7 

6.8 

Front  Wheel  Bounce,  u'nlj£ 

6.7 

6.6 

Rear  Wheel  Bounce, 

11 

10 

results  can  be  accepted  with  confidence.  This  verification  has  been  carried 
out  using  data  obtained  with  the  M-809  truck  during  field  trials  at  the  Aber- 
deen Proving  Grounds  (reference  3).  The  unsprung  mass  (rear  wheels) 
acceleration  has  been  chosen  for  comparison,  since  it  is  directly  related  to 
the  wheel  motions  from  which  the  shaker  signals  are  derived. 

The  result*  showing  the  PSD  of  this  parameter  for  the  point  contact, 
fixed  footprint  and  adaptive  footprint  tire  models  are  shown  in  Figure  11 
along  with  the  experimentally  determined  PSD.  As  can  be  seen,  the  point 
contact  model  significantly  overestimates  the  acceleration  over  the  entire 
frequency  range.  The  adaptive  footprint  model,  however,  shows  much 
better  agreement,  and  in  the  middle  and  high  frequency  bands  - of  particular 
importance  to  fatigue  simulation  in  the  shake  tests  - the  simulation  data 
closely  follow  the  experimental  results.  The  fixed  footprint  model  also 
shows  good  agreement  with  field  data,  although  at  the  natural  frequency 
peak,  it  somewhat  overestimates  the  force.  The  rigid  tread  band  model, 
though  providing  better  predictions  than  the  point  contact  model,  tends  to 
overestimate  the  force  in  the  middle  and  upper  frequency  range,  as  shown 
iu  Appendix  C. 
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Figure  11,  Comparison  of  Simulated  Response  with  Test  Data 
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These  reiults  confirm  the  initial  conclusion!  of  the  first  phase  of  this 
effort,  namely,  that  the  adaptive  footprint  tire  model  has  superior  capabili- 
ties for  predicting  the  transmitted  tire  forces,  and  can  be  used  with  advan- 
tage for  vehicle  simulation  and  shake  testing.  The  results  also  confirm  that 
the  fixed  footprint  tire  model  - though  predicting  slightly  higher  forces  in  the 
middle  frequency  band  - provides  a viable  alternative  for  simulations  in  which 
computation  requirements  are  to  be  minimised.  * 


♦As  explained  in  Appendix  A,  the  fixed  footprint  model  (with  linear  stiffness 
and  damping)  is  equivalent  to  a point  contact  model  moving  over  a filtered 
profile.  Profile  filtering  can  thus  be  carried  out  off-line,  and  the  dynamic 
simulations  can  then  proceed  using  the  simple  relationships  of  the  point 
contact  model. 
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3.  CONTROL  SYSTEM 


The  recorded  simulation  control  echeme  has  been  deacribed  in  terma 
of  two  atopa.  The  firat  atep,  obtaining  axle  displacement  records  from  a 
vehicle  simulation  program,  was  deacribed  in  the  laat  chapter.  Thia  chapter 
describes  how  the  axle  displacement  records  are  used  to  produce  input 
aignals  for  the  vehicle  shaker. 

There  are  several  subtaaka  involved  in  performing  the  second  atep. 

The  tasks  of  reading  the  axle  displacement  records,  storing  them  in  an 
easily  retrievable  form  and  supplying  them  to  the  shaker  at  a correct  rate 
demand  a sophisticated  control  system.  The  control  system  used  in  this 
effort  is  described  in  the  following  pages.  Detailed  procedures  for  operating 
the  control  system,  however,  are  not  described  in  this  rep  rt  since  they  are 
best  discussed  in  the  User's  Manual  (refersnce  4). 

3. 1 System  Requirements 

To  be  suitable  for  the  application,  the  control  system  must  meet 
several  requirements.  The  controller  must  be  able  to  synchronise  the 
shaker  input  signals  to  real-time,  otherwise , the  vehicle  will  not  experience 
the  inputs  in  the  same  time  frame  as  it  does  in  the  field  use.  To  be  able  to 
supply  the  data  in  real-time,  the  axle  displacement  records  must  have  a time 
scale  stored  along  with  the  displacement  values  and  a clock  should  interface 
with  the  controller  to  ensure  that  the  displacement  values  are  transferred  at 
the  stored  time  values.  For  this  application,  a high  frequency  clock  should 
be  used  so  that  the  accuracy  of  the  stored  time  data  is  preserved. 

The  data  transfer  rate  requirement  is  closely  coupled  to  the  previous 
requirement,  because  the  real-time  capability  can  be  achieved  only  if  the 
transfer  rate  is  adequate.  Data  transfer  rate  depends  on  tlms  step  between 
the  two  consecutive  displacement  valuea,  which  in  turn  is  derived  from  the 
maximum  frequency  content  the  signal  should  have.  The  time  step,  which 
is  fixed  while  running  the  vehicle  simulation  to  create  the  axle  displacement 
records,  should  be  such  that  the  displacement  signal  up  to  100  Ha  can  be 
reconstructed  from  the  digital  records.  This  requires  that  a time  step  of 
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about  5 milliseconds  should  be  used*.  The  displacement  data  at  this  time 
step  are  stored  in  16  bit  format  and  supplied  to  three  channels  (front,  middle 
and  rear  axle),  giving  a system  requirement  for  a minimum  data  transfer 
rate  of  about  9600  bits  per  second  (Baud).  If  the  upper  limit  on  the  frequency 
content  of  the  input  signal  is  increased  or  if  the  data  is  stored  in  32  bit 
floating  point  format  instead  of  16  bit  format,  the  transfer  rate  will  be  higher 
than  9600  bits  per  second. 

The  control  system  should  have  a capability  of  repeating  the  axle  dis- 
placement records  during  the  test.  Computing  and  storage  requirements 
generally  limit  the  record  length  to  less  than  a minute**  (it  is  18  seconds 
for  the  example  used  in  the  report).  At  the  end  of  the  record,  the  system 
should  be  able  to  start  from  the  beginning  without  a significant  loss  of  time. 

The  next  requirement  to  be  mat  by  the  control  system  is  adequate  data 
storage.  For  the  18  second  record  used  in  the  illustrative  example,  a 150 
thousand,  16  hit  work,  storage  medium  is  needed  to  include  axle  displacement 
information  along  with  the  corresponding  time  reference. 

The  final  requirement  for  the  control  system  is  a provision  for  con- 
verting the  stored  digital  data  to  analog  input  signals  for  the  shaker.  Digital 
to  analog  converters  with  bipolar  (signed)  outputs  and  at  least  11  bit  (0.  05 
percent)  accuracy  would  be  adequate.  The  digital  to  analog  conversion  rate 
should  also  be  high  enough  to  meet  the  required  data  transmission  rate. 

The  requirements  outlined  above  eliminate  from  consideration  the  use 
of  magnetic  tape  or  core  memory  of  a digital  computer,  because  the  magnetic 
tape  cannot  be  rewound  fast  enough  to  meet  the  requirement  for  repeating  the 
record,  and  the  core  memory  of  a computer  generally  does  not  have  a large 
enough  capacity  to  store  the  data. 


*The  smallest  frequency  that  can  be  discriminated  from  a digital  record 
with  time  step  T is  about  1/2T  Hz. 

**The  length  of  the  record  should  be  such  that  the  spectral  content  of  the 
profile  is  preserved.  See  reference  2 for  details. 
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• The  requirements,  then,  are  beat  met  by  a digital  or  hybrid  computer 

with  a magnetic  dlac  atorage  unit.  The  data  handling  and  timing  capabilitiea 
cf  a digital  computer  are  ideal  for  the  atorage  and  tranemiaaion  of  a large 
volume  of  data  aa  required  by  thia  ayatem.  A minicomputer,  with  adequate 
core  memory  to  run  the  control  software,  interfaced  with  a diec  atorage  unit, 
a digital  clock,  and  a bank  of  digital  to  analog  convertera  conatitutea  the  con- 
troller,  the  actual  hardware  of  which  ia  deacribed  in  the  next  section, 

3.2  System  Description 

The  control  system  selected  for  implementing  the  recorded  simulation 
control  scheme  la  an  EA1  (Electronic  Associates,  Inc. ) 590  hybrid  computer 
system  shown  schematically  in  Figure  12. 

Aa  ahown  in  the  figure,  the  central  part  of  the  system  ia  a Pacer  100 
minicomputer  with  a 16K,  16  bit  word,  core  memory.  The  digital  computer 
is  connected  through  a hybrid  interface  to  an  EA1  580  analog  computer  which 
incorporates  the  digital  to  analog  convertera  (D/A),  analog  to  digital  con- 
vertera (A/D),  sense  and  control  linea,  and  a 1 MHk  digital  clock  timer.  The 
digital  computer  ia  alao  connected  to  a 2. 2 million  word  (16  bit)  disc  atorage 
unit  through  a direct  memory  access  channel  (DMAC).  Input-output  to  the 
control  system  is  achieved  through  a high  speed  paper  tape  reader /punch,  a 
teletype  and  a line  printer.  Specifications  of  these  hardware  components  are 
summarised  in  Table  7. 

The  above  system  meets  the  requirements  outlined  in  the  previous 
section  quite  well.  The  clock  has  an  effective  frequency  of  42  KHa  (after 
taking  into  account  24  jiaec  required  by  the  read  and  set  subroutine!  con- 
trolling the  clock),  which  ia  adequate  for  synchronising  time  of  the  input 
signal  to  about  0.  5 percent  of  the  corresponding  value.  The  D/A  converters 
have  a setup  and  transfer  time  of  252  ms«c  for  the  three  channels,  thersfors, 
the  maximum  transmission  rate  is  close  to  4 KHa,  whioh  is  about  ten  times 
faster  than  required.  The  reset  time  for  the  system  after  reaching  the  end 
of  the  record  ia  about  0.  25  second,  which  is  considered  adequate.  Data 
storage  capacity  of  the  disc  ia  2,  2 million  words  as  compared  to  the  required 
capacity  of  150  thousand  words. 
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Pacer  100 


Figure  12.  Schematic  Diagram  of  the  Control  System  Hardware 


Table  7.  Control  System  Hardware 


PACER  100  Computing  System 

• 16,  384  1 6-bit  words  of  core  memory 

• Hardware  multiply  divide 

• Priority  interrupt  system 

e Memory  protect  system 

e Powor  fail/auto  restart  system 

e Five  universal  controller  positions 

e 1 Microsecond  memory  cycle 

Direct  Memory  Access  Channel  (DMAC) 

e Provides  for  I/O  concurrent  with  computation  on  cycle 
stoaling  basis 

e 1 Million  words /second 

Interval  Timer/ Real  Time  Clock 

e Provides  a high  resolution  16  bit  real-time  clock 

e Timer  update  signal  source  for  four  PACER  interval 
timers 

Disc  Cartridge  Controller  and  One  Dual  Disk  Drive  on  DMAC 

e Provides  1,  1 million  words  of  storage  on  one  fixed  platter 

• Provides  1,  1 million  words  of  storage  on  one  removable 
platter 

Hybrid  Device  Interface  693 

e Parallel  processor  control 

e Logical  control,  sensing  and  interrupts 
e Coefficient  device  setting 

e Data  monitoring /display 

e Eight  channels  D/A  multiplication  (15  bit) 
e IS  Channels  A/D  conversion  (13  bit) 

Input/Output  Devices 

• High  speed  paper  tape  reader,  300  characters/ sec. 
e High  speed  paper  tape  punch,  120  characters/sec. 

• Teletypewriter  station  (console  device),  ASR  33 

• Line  printer,  80  column 

• CPU  control  front  panel 


Table  7.  Control  System  Hardware  (Cont) 


Software  Systems 

• MHDOS,  moving  head  disc  operating  system 
e FORTRAN  IV  compiler 

e Assembler 

e BTE,  basic  text  editor 

e CIG,  core  image  generator 

e HOI,  hybrid  operations  interpreter 

e COP,  control  options  processor 

e FIU,  file  interactive  utility 

e RTL,  run  time  library 

e HLR,  hybrid  linkage  routines 

e Applications  software  packages 


Thus  the  control  system  meets  all  the  requirements  for  the  recorded 
control  scheme  application.  The  data  transfer  rate  and  the  storage  capacity 
are  significantly  more  than  those  presently  required.  But,  then,  room  is 
available  for  expansion  if  higher  frequency  input  signals  or  longer  record 
lengths  are  required, 

3.  3 System  Operation 

The  control  system  reads  and  stores  the  axle  displacement  records 
generated  by  the  vehicle  simulation  program.  The  records  are  transferred 
to  a buffer  and  converted  to  time  synchronised  analog  shaker  input  signals 
during  a vehicle  test.  This  simple  description  of  the  system  operation  is 
illustrated  by  Figure  13.  The  operation  described  in  the  figure  requires 
several  steps  by  the  user,  such  as  loading  the  program  into  memory,  open- 
ing and  positioning  the  data  file,  and  initiating  execution.  A complete  descrip- 
tion of  the  operation  ic  presented  in  the  User's  Manual  (reference  4)  and 
therefore  not  repented  here.  An  outline  of  the  operation  is,  however,  shown 
in  Figure  14  and  described  below. 
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Figure  13.  Control  System  Operation  Diagram 
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Figure  14,  Control  Program  Flowchart 
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The  operation  is  initiated  by  transferring  the  axle  displacement  data 
from  the  paper  tape  in  blocks  of  ten  time  steps  to  88  word  records  on  disc. 
This  format  of  storing  data  results  in  the  highest  rate  of  data  transfer* 

The  execution  begins  by  setting  the  control  pointers,  initialising  pro- 
gram parameters  and  activating  the  hybrid  interface*  The  clock  is  reset  and 
a data  block  containing  ten  time  steps  is  transferred  from  the  disc  to  a buffer 
in  cere  memory*  The  clock  time  is  then  compared  with  the  stored  time 
value.  If  the  clock  time  is  less  than  the  stored  time,  the  system  wait*  until 
it  becomes  equal  and  then  transfers  the  data  in  the  D/A  buffer  to  D/A  output. 
At  the  end  of  each  block  of  ten  time  steps,  more  data  is  retrieved  from  the 
disc  to  reside  in  the  buffer. 

At  the  end  of  each  block,  the  system  checks  for  the  end  of  the  axle 
displacement  record.  If  end  of  the  record  is  reached,  the  disc  is  rewound, 
the  clock  register  is  reset  and  the  process  is  repeated. 

3.4  System  Verification 

Performance  of  the  control  system  in  generating  shaker  input  signals 
from  axle  displacement  records  is  demonstrated  by  Figures  15,  16,  and  17. 
In  these  figures,  the  plots  of  the  three  axle  displacement  records  obtained 
from  the  vehicle  simulation  are  compared  with  the  outputs  of  the  respective 
D/A  channels.  For  a conversion  factor  of  10  volts  equals  1 foot,  the  plots 
are  identical  in  each  case,  thereby  verifying  the  system  performance. 


FP-ONT  AXLE  DISPLACEMENT  (ft. 


DISPLACEMENT  (ft. 


4.  VERIFICATION 


One  of  the  dominant  considerations  in  evaluating  merits  of  a method  for 
testing  vehicles  in  laboratory  is  correlation  of  the  inputs  to  the  vehicle  gen- 
erated by  the  test  method  with  those  experienced  in  field  environment. 
Translated  to  the  present  application,  this  means  that  the  axle  displacement 
inputs  produced  by  the  recorded  simulation  control  scheme  should  match  the 
axle  displacements  experienced  by  the  vehicle  for  similar  conditions  in  the 
field  environment. 

In  case  of  random  signals,  as  those  being  considered;  correlation  is 
best  expressed  by  similarity  of  the  spectral  densities.  Earlier  in  the  report 
the  spectral  densities  of  the  rear  axle  acceleration,  generated  by  the  vehicle 
simulation  program  employing  various  tire  models,  were  compared  with  the 
corresponding  data  obtained  in  the  field  test.  The  axle  displacement  spectral 
densities  for  the  same  test  are  now  compared  to  determine  if  the  vehicle 
experiences  similar  inputs  during  both  the  shake  test  and  field  test. 

In  monitoring  the  vehicle  response  during  the  field  test,  accelerometers 
were  used  to  measure  axle  acceleration  spectral  densities.  These  data  have 
to  be  converted  to  axle  displacement  spectral  densities  before  they  can  be 
compared  with  the  corresponding  input  signals  generated  by  the  control  sys- 
tem. Also,  the  measured  spectral  densities  have  to  be  plotted  against  the 
reduced  frequency  (£1),  and  not  frequency  in  Herts,  because  that  is  the  scale 
used  in  the  analytical  results,  These  transformations  are  achieved  through 
the  following  equations, 

S U>  „ s <*>  «2 

accln  in  accln  in 

ft/sec^  g 


s (0)  a S <*)  X- 

accln  accln  , - 
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S (0) 

; (0)  B accln 

aispln  b)  4 


where 

f = frequency  in  Hsb 

u » frequency  in  rad/ sec 

= reduced  frequency  in  rad/ft 

V * vehicle  velocity  in  ft/«ec 

g = acceleration  due  to  gravity 

S = spectral  density 

The  spectral  densities  of  the  front,  middle  and  rear  axle 
thus  obtained  from  the  field  data  are  then  plotted  along  with  the  corresponding 
spectral  densities  of  the  axle  displacement  records  produced  by  the  recorded 
simulation  control  scheme,  as  shown  in  Figures  18,  19  and  20.  The  figures 
show  that  the  spectral  densities  of  the  axle  displacements  obtained  from  the 
field  tests  are  in  good  agreement  with  those  generated  by  the  control  scheme. 
Just  as  discussed  in  Section  2,  employing  adaptive  footprint  model  and  fixed 
footprint  model  generally  result  in  more  accurato  representation  of  the  field 
environment  than  that  obtained  using  point  contact  or  the  rigid  tread  band 
model. 

For  all  three  axle  displacements,  the  agreement  amongst  the  models 
and  the  field  data  is  poorest  in  the  low  frequency  range.  In  that  range,  tho 
spectral  densities  are  calculated  based  on  only  a small  number  of  complete 
cycles.  Therefore,  there  is  a scatter  in  the  model  predictions.  If  the  profile 
length  is  increased,  the  scatter  will  reduce  and  the  models  will  agree  better 
with  the  test  data. 

For  the  front  axle  displacement  results,  the  adaptive  tire  model  is  able 
to  duplicate  the  field  test  data  quite  accurately  up  to  about  20  Hz.  At  higher 
frequencies  the  field  data  has  higher  power  content  than  that  predicted 
analytically.  This  is  partially  due  to  noise  in  the  signal  recording  system. 
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In  Figure  1%  the  advantages  of  employing  either  the  fixed  footprint 
model  or  the  adaptive  footprint  model  become  more  apparent.  Using  these 
models  produces  axle  displacement  signals  which  duplicate  the  field  axle 
displacements  quite  accurately,  whereas  using  the  point  contact  model  over- 
estimates the  spectral  density  by  a factor  of  10^  in  the  10-30  Hx  range. 
Similar  conclusions  can  be  drawn  from  the  results  shown  in  Figure  20,  (A 
spike  in  the  field  data  signal  in  this  case  at  about  60  Ha  may  be  due  to  the 
. noise  in  the  spectral  density  generating  system. ) 

A generally  good  agreement  between  the  analytical  and  the  field  test 
results,  shown  in  the  figures,  proves  that  the  scheme  is  able  to  provide  the 
vehicle  with  realistic  inputs.  Thus,  the  scheme  meets  its  development 
objective. 
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Figure  19.  Simulated  Response  - Middle  Wheel  Displacement 
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REAR  AXLE  D I SPIN 


5.  CONCLUSIONS 


The  recorded  simulation  control  scheme  is  developed  to  provide  a 
method  of  testing  vehicles  in  laboratory  without  requiring  field  data.  For  a 
particular  vehicle,  the  scheme  involves  simulating  vehicle  operation  over  a 
terrain  of  interest  at  a given  speed.  The  simulation  generates  records  of 
aide  displacements  which  can  be  played  back  to  provide  the  vehicle  shakers 
with  realistic  inputs.  Of  the  two  steps  involved  in  the  recorded  simulation 
control  scheme,  the  first  step,  i.  e.  producing  axle  displacement  records, 
uses  a terrain-tire-vehicle  computer  program.  The  second  step  is  per- 
formed by  a control  system  which  stores  the  axle  displacement  records  and 
converts  them  to  time- synchronized  analog  signals  during  a shake  test. 

The  vehicle  simulation  program  incorporates  four  different  models  of 
tiros.  Applicability  of  these  tire  models  is  explored  by  comparing  the  simu- 
lation results  with  those  obtained  during  field  testing  of  an  M-B09  cargo 
truck.  The  comparison  shows  that  the  adaptive  footprint  model  provides 
superior  tire  force' predictions,  due  to  its  ability  to  envelop  ground  irregu- 
larities through  local  footprint  deformations.  Vehicle  simulations  using  the 
adaptive  footprint  model  are  well  able  to  predict  the  principal  natural  fre- 
quencies of  vehicle  vibration,  and  the  simulated  response  agrees  closely  with 
field  trial  data,  particularly  in  the  middle  and  upper  frequency  range.  Since 
many  shake  test  failures  are  caused  by  high-frequency  stress  reversal 
(fatiguo),  the  adaptive  footprint  model  in  particularly  well- suited  for  simu- 
lating fatigue  failures  in  the  laboratory. 

The  fixed  footprint  model,  while  not  as  advanced  as  the  adaptive  foot- 
print model,  also  provides  improved  force  predictions.  Since  this  model 
can  be  simulated  by  an  equivalent  point  contact  model  in  series  with  a filter 
(see  Appendix  A),  it  will  be  a viable  alternative  for  simulations  in  which  the 
computational  requirements  are  to  be  minimized.  The  simpler  tire  models 
(i. e.,  the  point  contact  and  rigid  tread  band  models)  generally  overestimate 
the  tire  forces,  particularly  at  high  frequencies  (1,  e, , typically  in  the  5-50 
Hz  region),  and  their  use  will  be  governed  by  this  limitation. 
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The  control  system  perform*  It*  function  adequately , The  ayetem  1* 
able  to  atore  the  axle  dieplacement  record*  obtained  from  the  vehicle  simu- 
lation program  on  a magnetic  disc  and,  at  the  time  of  shake  teat,  produce* 
continuous  analog  input  signals  by  repeating  the  limited  length  records. 
These  analog  signals  are  synchronised  perfectly  to  the  real  time  and  there 
is  no  detectable  difference  between  the  analog  output  and  the  digital  input  to 
the  control  system. 

The  axle  displacement  signals  obtained  are  compared  to  the 
corresponding  signals  recorded  in  a field  test.  Both  signals  show  similar 
spectral  densities,  thereby  proving  that  the  recorded  simulation  scheme  is 
able  to  provide  adequate  control  to  the  vehicle  shaker, 


6.  RECOMMENDATIONS 


The  recorded  simulation  control  scheme  described  in  this  report  can 
be  used  in  many  ways  to  improve  performance  of  vehicles.  The  scheme  can 
be  employed,  for  example,  to  carry  out  a study  of  fatigue  failures  of  vehicle 
components  by  shake  testing  the  vehicle  for  extensive  periods  of  time.  The 
tests  can  be  repeated  for  different  terrains  and  speeds  and,  from  the  test 
results,  modifications  to  prevent  such  failures  can  be  suggested. 

A study  to  select  optimum  tireB  for  a particular  application  can  be 
performed  using  the  control  scheme.  In  such  a study,  data  for  candidate 
tires  can  be  obtained  from  manufacturers  and  performance  of  the  vehicle 
incorporating  those  tires  can  be  determined  using  the  vehicle  simulation 
program  and  the  shake  test  facility.  A comparison  criteria,  such  as  the 
maximum  speod  at  which  the  driver  motion  becomes  unacceptably  large, 
can  evaluate  performance  of  various  tires  and  the  optimum  tire  can  be 
selected. 

The  recorded  simulation  control  scheme  can  also  be  used  to  carry  out 
a parametric  analysis  of  vehicles.  Effects  of  changing  various  parameters, 
such  as  the  suspension  spring  rates,  damping  coefficients,  center  of  gravity 
location,  etc,  , on  vehicle  performance  can  be  predicted  by  the  computer 
simulation  program  and  verified  by  the  shake  test  facility.  Such  a parametric 
analysis  can  be  helpful  in  designing  new  vehicles  or  modifying  the  existing 
ones. 


The  recorded  simulation  scheme  itself  can  be  modified  to  expand  its 
capabilities.  For  example,  the  core  memory  incorporated  in  the  present 
control  system,  which  is  sufficient  to  perform  record  storing  and  input 
signal  supplying  steps  of  the  scheme,  can  be  increased,  so  that  the  vehicle 
simulation  program  can  be  executed  on  the  control  system  and  the  inter- 
mediate step  of  preparing  axle  displacement  records  on  paper-tapes  (or 
some  such  medium)  can  be  eliminated.  The  control  system  can  also  be 
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modified  to  record  outputs  of  various  transducers  monitoring  performance 
of  the  vehicle.  These  records  can  then  be  processed  to  obtain  statistical 
information,  such  as  spectral  densities  and  cross -cor  relation  of  the  variables 
describing  vehicle  performance. 

The  vehicle  simulation  program  can  be  expanded  to  include  roll  motion 
of  the  vehicle.  Also,  the  fore-and-aft  forces  generated  by  the  simulation  can 
be  used  in  vohicle  testing  if  the  vehicle  shaker  is  suitably  modified. 

A final  suggestion  for  future  work  is  to  develop  synthesized  simulation 
control  (described  in  Chapter  1)  because  this  scheme  eliminates  the  task  of 
storing  the  axle  displacement  records.  Also,  by  using  random  signal 
generators,  biases  which  are  present  in  the  stored  records  are  removed. 
However,  developing  such  a scheme  may  need  significant  software  and/or 
hardware  development.  Meanwhile,  the  recorded  simulation  control 
scheme  can  provide  an  improved  method  for  laboratory  testing  of  vehicles. 
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APPENDIX  A 


THE  TIRE  MODELS 


Details  of  the  four  tire  models  are  presented  in  the  following: 
A-l.  Point  Contact  Tire  Model 


The  schematic  diagram  of  the  point  contact  tire  model  is  shown  in 

Figure  A-l.  The  terrain  elevation  coordinate  yQ  is  measured  positive 

upwards  from  the  mean  line  and  is  a function  of  the  distance  x.  The  wheel 

center  motion  coordinate  U measured  positive  upwards  from  the  initial 

equilibrium  position:  i.  e. , with  the  tire  statically  loaded  under  the  vehicle 

weight  and  yQ  = 0.  Fy  and  F^  are  the  vertical  and  horizontal  components 

of  footprint  force.  F and  F are  the  tire  forces  transmitted  to  the  vehicle 

y x 

suspension.  The  principal  assumptions  of  the  analysis  are: 

1.  The  tire  mass  is  concentrated  at  the  wheel  center. 

Z.  Terrain  contact  occurs  at  a single  point  vertically  below  the 
wheel  center. 

3,  The  total  footprint  force  acts  normal  to  the  terrain  surface  at 
the  contact  point, 

4,  When  in  contact,  forces  between  the  follower  and  the  ground 
must  be  compressive.  When  the  contact  force  becomes  zero, 
the  follower  leaves  the  ground  (wheel  hop). 
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Direction  of  Motion 


Resultant  Tire  Force 
Transmitted  to  Suspension 
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Figure  A-l.  Schematic  Diagram  of  Point  Contact  Tire  Model 
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'*8*1 


5.  The  tire  force. deflection  characteristics  are  represented, 
in  general,  by  nonlinear  stiffness  k and  damping  b.  In 
general,  k and  b can  be  functions  of  the  tire  deflection 
and  deflection  rate. 


For  arbitrary  vertical  displacements  and  velocities  y^,  y^,  , 

yj  the  total  spring-damper  force  is 


' s 

r = / 


V*i 

kdy  + J*  bdy 


(A-l) 


where  y ^ is  the  static  deflection  of  the  tire  under  the  equilibrium  load 

W,  and  y is  the  time  rate  of  change  of  profile  elevation  seen  by  the 

follower  due  to  tire  motion,  y , and  y are  found  from 

'st  ‘o 


rBt 

i 


kdy  = W 


(A-2) 


Y0  = V(dyQ/dx) 


(A- 3) 


where  V is  the  forward  velocity  (dx/dt)  and  (dyQ/dx)  is  the  slope  of  the 
ground  at  the  contact  point. 


The  vertical  footprint  force  Fy  will  be  equal  to  F when  the 
follower  is  in  ground  contact  (F  > 0)  or  zero  when  the  follower  is  out 
■ of  contact  (F  0) 


F. 


} F > 0 


F < 0 


(A-4) 
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Since  the  resultant  footprint  force  acts  normal  to  the  terrain,  the  fore- 
and-aft  force  component  is  related  to  the  vertical  component  as 
follows. 


Fh/Fv 


dy  /dx 
o 


(A- 5) 


Equations  (A-l)  to  (A-S)  determine  the  footprint  forces  F and  F,  . The 
tire  forces  transmitted  to  the  vehicle  suspension,  F and  are  deter- 
mined from  the  equations  below. 


Fx  * Fh  (A- 7) 

A- 2.  Rigid  Tread  Band  Model 

The  schematic  diagram  of  the  rigid  tread  band  model  is  shown 
in  Figure  A-2.  Comparison  with  the  point  contact  model  of  Figure  A-l 
shows  that  it  is  equivalent  to  a point  contact  model  traversing  a modified 
terrain  profile  yQ(x)  obtained  from  the  tread  band  center  motion.  There- 
fore, the  analysis  and  equations  for  the  point  contact  model  shown  in 
the  previous  section  remain  unchanged,*  and  it  is  only  necessary  to  find 
the  modified  profile  yQ(x)  in  terms  of  the  original  terrain  profile  yQ(x) 
and  the  tread  band  radius  r. 


Because  of  the  presence  of  the  circular  tread  band,  the  terrain 
contact  point  will  in  general  be  shifted  a distance  x fore  or  aft  of  the 
wheel  center  as  illustrated  in  Figure  A-2.  From  geometry,  the  tread 
band  center  height  y^  at  any  location  x,  is  given  by 


y0(*) 


yQ(x  + x)  + 


(A-B) 


The  criterion  for  wheel  hop  also  remains  unchanged,  i.  e.  , when  the 
contact  force  goes  to  zero,  the  wheel  leaves  the  ground. 
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Direction  of  Motion 


Resultant  Tire  Force 
Transmitted  to  Suspension 


The  criterion  for  determining  x is  that,  at  the  contact  point, 
the  elope  of  the  tread  band  and  the  terrain  profile  must  he  equal.  The 
terrain  profile  slope  at  the  contact  point  is  d(y^(x  -I-  x)J  /dx.  The  slope 
of  the  tread  band  at  this  point  is  - d ( Vr^  - x^ ) /dx.  Equating  these 
two  slopes,  the  condition  that  must  be  satisfied  at  the  contact  point  is 


(A-9) 


Solution  of  the  above  equation  determines  the  contact  offset  distance  x 
in  terms  of  the  tire  location  x and  radius  r.  In  general,  the  roots  of 
Equation  (A-9)  will  give  two  values  of  x corresponding  to  diametrically 
opposite  points  having  equal  slope,  as  shown  in  Figure  A- 3.  Physically, 
contact  can  only  occur  in  the  lower  half  of  the  tread  band,  and  for  this 
point  dG/dx  will  be  negative.  Solutions  of  Equation  (A-9)  can  also  give 
multiple  pairs  of  x values  implying  that  the  slope  condition  can  be  satis- 
fied at  several  points  on  the  lower  half  of  the  tread  band.  In  such 
cases,  the  highest  value  of  yQ(x)  found  from  Equation  (A-8)  will  repre- 
sent the  physically  realizable  condition.  The  other  values  can  only  be 
satisfied  through  tread  band  deformations  as  illustrated  in  Figure  A- 3. 

A section  of  modified  profile  obtained  as  described  above  through  wheel 
center  filtering  is  included  in  Figure  7.  The  original  profile  from  which 
the  modified  profile  was  derived  is  also  included  in  this  figure. 

Thus  Equations  (A-8)  and  (A-9)  along  with  the  criteria  for  selecting 
the  physically  admissible  root  allow  determination  of  the  tread  band  center 
motion  yo(x),  which  along  with  Equations  (A-l)  through  (A -7) 
constitute  the  rigid  tread  band  analytical  model. 


A- 3.  Fixed  Footprint  Model 


The  schematic  diagram  for  the  fixed  footprint  model  is  shown 
in  Figure  A-4.  In  this  model,  terrain  contact  occurs  ever  a finite 
footprint  length  L,  and  the  tire  stiffness  and  damping  is  uniformly 
distributed  over  the  contact  length.  The  principal  assumptions  of  the 
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(a)  Physical  Representation  of  Roots  of  Equation  (A- 7) 


Figure  A-3.  Determination  of  Admissible  Terrain  Contact  Point 


Direction  of  Motion 


I**—  Footprint  Length,— +-| 
L 


Figure  A-4.  Schematic  Diagram  of  Fixed  Footprint  Tire  Model 


analysis  are  summarized  below. 


1.  The  tire  mass  is  concentrated  at  the  wheel  center. 

2.  Terrain  contact  occurs  through  a footprint  of  fixed  size, 
independent  of  tire  deflection. 

3.  The  footprint  center  is  constrained  to  lie  vertically 
below  the  wheel  center. 

4.  The  footprint  force  components  act  normal  to  the  local 
terrain  aurfaco. 


5.  When  in  contact,  the  resultant  footprint  pressure  force 
must  be  comprosiive,  and  all  parts  of  the  footprint 
must  contact  the  ground.  When  this  force  becomes 
aero,  the  complete  footprint  leaves  the  ground. 

6.  The  tire  fore -deflection  characteristics  are  represented, 
in  general,  by  nonlinoar  stiffness  and  damping  uniformly 
distributed  over  the  footprint. 


For  arbitrary  vertical  displacements  and  velocities  y (x),  y,, 

* • O 1 

yQ(x).  y^  the  total  force  in  the  distributed  spring -damper  olemonts  is 

+ ii/2  y,t  + y0®-vj 

F a f I k' dy  dx 

?L/2  Jo 


+ hU 

+ f f b1  dy  dx 


(A.  10) 


where  y it  the  (uniform)  static  deflection  of  the  footprint  under  the 
equilibrium  load  W,  yQ(x)  is  the  profile  elevation  as  a function  of  the 
footprint  length  coordinate  x,  and  y0(x)  is  the  time  rate  of  change  of 
profile  elevation  at  any  location  x within  the  footprint  due  to  tire  forward 
motion,  y ^ and  yQ(x)  are  found  from 


y.t 

f k'Ldy*  W 


(A- 11) 


. _ d yQ<x) 

y (X)  = V 

° dx 


(A-12) 


where  V is  the  forward  velocity  and  dyQ(x)/dx  is  the  local  slope  of  the 
ground  at  the  footprint  location  x. 

As  discussed  in  Section  A-l,  for  the  point  contact  model,  the 
footprint  force  Fy  is  given  by 


F * F 
v 


F » 0 

V 


F > 0 


F < 0 


(A- 13 ) 


d y (x) 

rh/rv  ■ -tS— 


(A- 14) 


Equations  (A«*10)to  (A-14)  determine  the  footprint  forces  Fy  and  for 
the  fixed  footprint  tire  model.  The  tire  forces  transmitted  to  the 
vehicle,  F and  F^,  are  determined  from  the  basic  tire  model  equations, 
Equations  (A- 6)  and  (A-?). 
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Fy  “ Fy  “ mVl 


F - F, 
x h 


For  the  special  ease  of  linear  stiffness  and  damping,  the  fixed 
footprint  tire  model  is  equivalent  to  a point  contact  model  having  the 
same  total  stiffness  and  damping  (i.  e.  , k = k'L,  b = b'L)  and  travelling 
, over  a modified  terrain  profile  such  that  each  point  on  the  modified 
profile  is  the  mean  value  of  the  original  profile  taken  over  the  footprint, 
i.  e. , 


+ L/2 


yc(x) 


* / 


yo(x  + x)  dx 


(A- 15) 


This  is  equivalent  to  reducing  the  amplitude  of  every  Fourier  component 
of  the  original  profile  by  a factor  of  Sin{TrL/A)/(irL/X)  where  A is  the 
wavelength  of  the  component  considered.  ' A section  of  filtered  profile 
obtained  as  described  above  5s  included  in  Figure  7.  The  original 
profile  from  which  this  was  derived  is  also  shown  in  the  figure. 

A-4.  Adaptive  Footprint  Model 

The  schematic  diagram  for  the  adaptive  footprint  tire  model  is 
shown  in  Figure  A-5.  In  this  model,  terrain  contact  occurs  over  a 
footprint  of  finite  length,  which  is  dependent  on  tire  deflection.  Tire 
stiffness  and  damping  is  uniformly  distributed  over  the  angular  segment 
that  includes  the  contact  zone.  The  principal  assumptions  of  the  analysis 
are  summarized  below. 

1.  The  tire  mass  is  concentrated  at  the  wheel  center. 


Every  tread  element  deforms  independently  of  its  neighbors 
and  contac  a the  terrain  in  the  footprint  zone.  The  tire 
retains  its  original  shape  outside  the  footprint. 


! 

9 

j 

i 

I 

3.  The  tire  force>deflection  character!  a tics  are  represented 
as  follows: 


(a)  By  a constant  inflation  pressure  p^  acting  over  the 
footprint  area. 

1 

i 

(b)  By  distributed  nonlinear  radial  stiffness  k"  and 

damping  b"  to  simulate  carcass  load  contributions. 

* ! 

As  in  the  other  models,  k"  and  b"can  be  functions 
of  the  tire  deflection  and  deflection  rate. 


1 

. _ . 1 

For  arbitrary  vertical  displacements  yQ(x),  y^»  yo(x),  Yj*  the 
force  components  d and  d due  to  deflection  of  a tread  element  j 

oriented  at  an  angle  6 to  the  vertical  (initial  distance  x from  the  wheel 
center)  is  given  by  4 


d F (0)  s d F cose  + 
v c 


Pi  B rdfi 


/TT 


d yQ(x) 

f — i. 

L dx 


d Fv(e)  = 0 


dF  > 0 
c 


d F < 0 

C *■“ 


j 


(A-  16) 


d Fh(e)  = d Fc  Sine  + 


d Fh(6)  = 0 


dFc>  0 


d F < 0 
c — 


(A-  17) 
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H 


where 


r6(9)  6(9) 

* \ j k"  dy  + C b"dy 


(A- 18) 


x = r Sin  9 


(A- 19) 


la  Equations  (A-l6)  and  (A- 17),  the  first  term  repreaente  the  carcaaa 
force  contribution,  and  the  second  gives  the  inflation  pressure  force. 
In  these  equations 

B - effective  footprint  width 


dyo(*) 


local  slope  of  terrain  profile  at 
contact  point 


radial  deflection  of  tread  element 
orionted  at  angle  0 from  vertical 

radial  velocity  of  tread  element 
due  to  tire  motion. 


The  radial  deflection  6(0)  is  the  sum  of  the  initial  deflection 
68t(0)  due  to  the  equilibrium  load  W,  and  the  deflection  due  to  terrain 
irregularity  yQ(x)  and  wheel  center  motion  y^ 


6(0)  = 6>fc(0)  + yQ(x+x)  - yj  /Cose 


'dy0(«4S)j  _ ^ | 


4(8)  = 


/Cos  6 


(A-  2 q 


(A-21 ) 
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The  initial  radial  deflection  6# fc( Q)  at  location  e under  the  static 
vehicle  weight  is  found  from  solution  of  the  following  three  simultaneous 
equations 


/ 


cos 0k”  dy  d0  + A 


W 


6.t<0) 


yBt  “ “ Co8  0 ) 

Cos  9 


(A-22) 


(A-23) 


(A-24) 


where  y#t  is  the  equilibrium  tire  deflection  at  the  footprint  center  (0  = 0), 
and  A is  the  effective  contact  area  at  equilibrium. 

The  effective  tire  width  B is  equal  to  the  effective  area  divided 
by  the  footprint  length, 

B = A/2r  Sin  9 (A-2!>) 

o 

* 

Finally,  the  total  vertical  and  horizontal  components  of  footprint 
force,  F and  F^,  can  be  obtained  by  integrating  Equations  (A-  16;  and 
(A-  1 7)  over  the  lower  half  of  the  tire,  where  terrain  contact  is  possible 

9=  +ir/2 

Fv  = / dFv(9)  (A“26) 

0=  - ir/2 

Z + n/2 

dFh(9)  (A- 2 7) 

- W 2 

The  forces  transmitted  to  the  vehicle  suspension,  F and  F , can  then 

y x 

be  found  from  Equations  (A-6)  and  (A-7)  as  before. 


APPENDIX  B 


THE  VEHICLE  MODEL 


The  vehicle  model  includes  all  the  vehicle  components  excluding  the 
tires.  The  tire  model  is  derived  separately  (see  Appendix  A)  and  the  general 
tire  model  variables  are  converted  to  the  vehicle  model  variables  as  shown 
in  Table  B-l.  The  vehicle  model  consists  of  the  following  elements:  the  hull, 
the  bogie,  the  front  suspension,  the  middle  suspension,  the  rear  suspension, 
the  front  unsprung  mass,  the  middle  unsprung  mass  and  the  rear  unsprung 
mass.  These  models  are  described  below. 

B-l  The  Vehicle  Body 

• 

The  vehicle  body  or  hull  is  shown  in  Figure  B-l.  It  is  acted  upon 
by  horizontal  and  vertical  forces  from  the  front  suspension,  bogie  con- 
vection and  frame  mounted  stops  as  shown.  The  forces  are  assumed  to 
act  in  the  horizontal  and  vertical  directions  in  an  inertial  frame  of 
reference,  and  not  in  the  vehicle  reference  frame. 

Two  variables,  incremental  vertical  displacement  of  the  hull  CG 
(Yq)#  measured  from  initial  equilibrium  condition,  and  the  hull  pitch 
angle  (0)  defines  the  position  of  the  hull  at  any  instant. 

The  equations  of  motion  about  the  CG  are  as  follows. 

Mh  Yo  = “ “h*  + 2 <Fy2  + Fy34  * Ffa  + Fbs>  <>> 

^0=2  Fy34  [L3  cos  © + ( S34  + GG)  sin  0] 

+ 2 Ffs  [(L3  - Sf8)  cos  G + GG  sin  © ] 

+ 2 Fbs  t(L3  + Sbs)  cos  © + GG  sin  ©] 

+ ?.  Fx34  [(S34  + GG)  coa  © - L3  sin  ©J 

• 2 Fy2  (L2  cos  0 - GG  sin  ©)  - Fy2  (L2  sin  © + GG  coa  9)  (2) 
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Hull  Weight  - M^g 

Front  Suspension  Forces  - F 

Bogie  Connection  Forces  - F 

Frame  Stop  Forces  - F F 

fs*  i 


Figure  B-l.  Freebody  Diagram  of  Hull 


The  factor  of  2 in  the  above  equations  appears  because  identical  forces 
act  on  the  right  and  left  side  of  the  vehicle,  and  the  forces  shown  in 
Figure  B-l  (except  the  weight)  are  for  one  side  only. 

B-2  Bogie 

The  bogie  is  pivoted  to  the  hull,  so  that  it  can  rotate  relative 
to  the  hull.  The  bogie 'hull  connection  is  considered  to  be  a frictionles 
pin -joint,  transmitting  only  forces  and  no  torques. 

The  bogie  mass  and  inertia  includes  those  of  the  leaf  springs 
and  pivot.  The  middle  and  rear  axle  masses  are  included  in  the 
unsprung  mass,  and  not  in  the  bogie  mass. 

The  freebody  diagram  (inertial  reference  frame)  of  the  bogie  is 
shown  in  Figure  B-2. 

The  variables  which  define  the  bogie  position  are  the  vertical 
displacement  of  the  bogie  CG  measured  from  the  initial  equilibrium 
condition,  (Y^),  and  the  bogie  pitch  angle  [(f)).  The  equations  of  motion 
for  vertical  motions  and  rotation  about  the  CG  are 

^ Yb  = - M,  g + 2 (Fgy3  + F#y4  - Fy34) 

% # - 2 (Fiy4  * cos  </)  + F-x4  a sin  0 + Fy34  ^ sin  0 - 
Fx34  Sfa  cos  0 ’ F.y3  a co>  0 ’ Fsx3  a *ln  $ 

From  a horizontal  force  balance, 

Fx34  5 Fsx3  + Fbx4 

Because  the  bogie  pivot  is  attached  to  the  hull,  the  bogie  CG 
« • • * • « 

acceleration  (Y^)  can  be  expressed  in  terms  of  Y , ® and  tp  . From 
Figure  B-3, 


I 

=1 


Direction  of  Motion 


Fsy3  ^8 


Rear  Suspension  Forces  - F p 
Middle  Suspension  Forces  - "p4  f "/4 
Bogie  Connection  Forces  - F p 8/3 
Bogie  Weight  - ^g  *34’  Y 34 


l 

l 

l 

f 

) 

I ; 

Figure  B-2.  Freebody  Diagram  of  Bogie 
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Figure  B-3.  Hull -Bogie  Linkage 


= Yq  + ■in  © (GG  + S^)  (1  - co»  ©)  + (1  - co»  <p) 

Differentiating Equation  (6)  twice 


(*>) 


Yb  = Yq  + L3  (-  ©2  sin  © + © cos  ©) 

+ (S^  + GG)  (®2  co«  © + © tin  ©) 

+ Sjj  (i2  cos  <p  + <p  Bin  <p)  (7) 

The  bogie  force  Fy34  can  be  determined  by  solving  Equation*  (1)  to  (4) 
along  with  Equation  (7)  simultaneously.  This  gives 


y34 


where 


/ 2XD 

XA 

2XG-  XC 

2XH.  XF 

• - XlJ 

te* 

“ “h  ‘ 

*b 

■ 

/ 2XG.  XB 

2XH-  XE 

2 . 

2 \ 

\ 

*h 

“b 

j 

XA 

* *t.  + Fb. 

) 

XB 

= 

COB  © + (S.j4 

+ GG)  sin  © 

xc 

it 

K 

• 

Kl3  - s(l) 

cos  © 4 GG 

■in  ©3 

+ Fb. 

Ul3  * 

co.  0 4 

• GG  sin  ©] 

“ Fx34 

[(S34  + GG) 

- L3 

sin  ©3  - F, 

yZ  <L2  CO“  0 

- GG  sin  ©)  - Fx2 

(L^  sin  6 + GG  co*  ©) 

(8) 


© 


ao 


g 


XD 

XE 

XF 


XG 

XH 

XI 


= F.y3  + Fsy4  “ *T 

= % <t> 

= Fsy4  a coa  + Fgx4  * ain  $ " Fx34  \ co*  0 

* FSy3  * COB  4>  ’ F8X3  6 8in  <*> 

= coa  0 + sin  © + GG  sin  © 

= ^ ain  0 

- - Li^  ©^  sin  Q + S,j4  ©^  coa  © + 0 ^ coa  0 

+ GG  coa  © 


B-3  The  Suspension 

There  are  two  pairs  of  leaf  springs  connecting  the  chassis  to  the 
axles.  The  front  pair  support  the  front  axle,  and  along  with  shock 
absorbers  form  the  front  suspension.  The  rear  leaf  springs  support  the 
middle  and  the  rear  axles  and  form  the  middle  and  the  rear  suspensions. 
There  are  no  shock  absorbers  in  the  middle  and  rear  suspensions.  The 
suspension  spring  and  damping  characteristics  are  shown  in  Figure  4. 


B-3. 1 The  Front  Suspension 


The  front  suspension,  shown  schematically  in  Figure  B-4, 
consists  of  the  leaf  spring,  modelled  as  a linear  spring  with  elastic 
stops  in  parallel  with  dry  friction  damping,  and  a shock  absorber.  The 
forces  on  the  suspension  are  as  follows! 


k2 


+ F 


b2f 


+ F,  , + F . , 

b2v  stop  2 


(9) 
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Direction  of  Motion  Fy2 


Fisturo  R-4.  Froebody  Diagram  of  Front  Sunpanaion 
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^k2”  spring  shii'fhess  force 


where 


where 


where 


6 


’■  K2  *2 

2 = Y.p2  + Vw2  ' Y.2 

- equilibrium  suspension  displacement 


- Mh8  L3 
sp2'  2 K2  (L2  + L3) 


Ys2  = Yo  " L2  sin  0 
Fb2£  -dry  friction  force 

*b»  5 Fr.f  K2  I Yw2  - Y.2I 

Fr«f  -dry  friction  factor 


(assumed  to  be  0,  05) 


Fb2v 

Fb2v 


*2  = Yw2  “ Ys2 

Ys2  S Yo  - L2  C0S  ® 
- shock  absorber  force 

= B2v  *2 


© 


B 


2vj 


2vr 


for  6 g 
for  6Z 


> 


< 


0 

0 


(10) 

(11) 

(12) 

(13) 

(14) 


(15) 

(16) 

(17) 
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atop  force 


Fstop  2 


F«top  2 


Y.2-JV°K2'  Yw2  * Y.2  >-  J, 


n2 


~ Rb2  < Yw2  “ Y.Z  < Jn2 


" Yw2  'Rb2>aV  Y»2  Yw2  ~ Rb2 


(18) 


It  should  be  noted  that  the  dry  friction  force  given  by  Equation  (14) 
is  approximate;  a more  advanced  model  requires  simultaneous  solution 
of  several  equations  and  has  not  been  included  in  this  initial  formulation. 

B-3.  2 The  Middle  and  Rear  Suspensions 

The  middle  and  rear  suspensions,  formed  by  the  bogie  leaf 
springs  are  identical.  Figure  B-5  shows  the  schematic  diagram  of  these 
suspensions. 


For  the  middle  suspension,  the  forces  are 


F.y3  8 Fk3  + Fb3f  + Fstop  3 


(19) 


Fk3  -spring  stiffness  force 


whore 


Fk3  = K3  63 


6_  « Y - + Y . - Y , 

3 sp3  w3  s3 


(20) 

(21) 


sp3 


Ysp3*  e qxsllib r ixim  suspension  displacement 

te-Tcrba  • ^ .a. 


(22) 
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Direction  of  Motion 


Figure 


B-5.  Freebody  Diagram  for  Middle  and  Rear  Suspensions 


Y.3 

= Y^  - a sin 

(23) 

II 

>* 

Y + L,  sin  © 
o 3 

+ (1  ■ cos  ©)  + 

®b  (1  - 

cos  <p) 

(24) 

Fb3f 

- dry  friction  force 

1 

. («,) 

Fb3f  8 Fr.f  K3  IYw3  " Yb3  1 'JJ’ 

(25) 

where 

*3  ■ Yw3  ' Y, 

» 3 

(26) 

Y.3  • Yb  ' » 

COS  <f)  <p 

(27) 

Y.  ■ Y + L- 
b o 3 

• 

cos  © © 

(28) 

+ sin 

© © + sin  0 0 

F.top  3 - ,t0P  £orce 

| 

< <Yw3  - Y.3  * 

j»3)  »K3 , yw3 

' Yb3  — 

Jn3 

Fatop  3 ° \ 

i 

L 

1 

' Sb3  < Yw3  ' 

■ Yb3  < 

Jn3 

(29) 

1 

^ - <Y.3  - Yw3 

- Rb3>  “K3>  Y.3-Yw3 

> Hb3 

Similarly  for  the  rear  suspe  iSion,  the  forces  are 


Fsy4 

S Fk4  + Fb4f  + Fstop  4 

(30) 

Fk4 

= K3  *4 

(31) 

11 

‘O 

YsP4  + Yw4  - Ys4 

(32) 

i- 

f 
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8p4 


. [ “h8  L2  . “b  « 1 1/K 

['4  (L2  + L3)  4 J 


Y*4  a Yb  + a «in  0 


Fb4f  “ Frsf  K3  IYw4  • Y.4>  34  / l64I 


S4  = Yw4  * Y.4 


Y«4  * Yb  + » “•  * <#. 


<Ywi  * Y.4  ' J“3>  ”K3‘  Yw4  - Y.4  2 *>5 


atop  4 


- “3  < Yw4  * Y.4  < J”3 


<Y.4  ' Yw4  ' Hb3>  nKV  Y.4  ' Yw4  > Eh 


s4  w4 


(33) 

(34) 

(35) 

(36) 

(37) 


(38) 


B-4  The  Unsprung  Maes 

The  front  unsprung  mass  (each  side)  consists  of  the  mass  of  half 
the  axle  and  one  front  tire,  while  the  middle  and  the  rear  unsprung  mass 
(each  side)  consist  of  the  mass  of  half  the  axle  and  two  tires  (because 
the  rear  axles  have  dual  tiree). 


B-4. 1 Front  Unsprung  Mass 


The  free  body  diagram  of  the  front  unsprung  mass  is 
shown  in  Figure  B-6(a).  The  equations  of  motion  in  the  vertical  and 
horizontal  directions  are 


^ Yw2  - - M,  g + Fv2  - F, 


y2 


Fx2  = FhZ 


(39) 


(40) 


Tirs  Tore** 


Suspension 

Forces 


FvZ*  Fv3'  Fv4 


FM*  Fh3*  Fh4 
Fx2*  Fix3*  Fsx4 


Fy2;  F.y3'  F.y4 


-J 


Unsprung  M u M g M g (*)  Front  Unsprung  Muss 

Weight  3B'  4* 


TU  {Fbu\  l^yS^syd* 


F*x3  ^Fsx4* 


j^wl  ^w4^ 


Fh3  <Fh4> 


Fv3  (Fv4' 


Mjg  (M^) 


Notei  Terms  in  { ) correspond’ 
to  reer  unsprung  mess* 


(b)  Middle  end  Resr  Unsprung  Mess 


Figure  B-6.  Freebody  Disgrsm  of  Unsprung  Mess 


Acceleration  terms  ere  absent  in  Equation  (40)  because  constant  forward 
velocity  has  been  assumed. 


i 


£-4.2  Middle  and  Rear  Unspruns  Mass 

The  free  body  diagram  of  the  middle  and 

rear  unsprung 

mass 

is  shown  in  Figure  B-6(b), 

s 

[ The  equations  for  the  middle  unsprung  mass  are: 

M3  *w3  * ’ M3  « + F,3  ' F.y3  ‘ Ff. 

(41) 

1 

1 

F.x3  * Fta3 

(42)  j 

1 

Similarly  the  equations  for  the  rear  unsprung  mass  are 

$ 

j 

| 

*„4  ■ - m4  « + - r.y4  - rb. 

(43)  | 

F«4  * FM 

(44) 

J 

Stop  forces  F^g  and  F^g  are  determined  as  follows: 

< 

1 

Cb  ■ Sb.  («-♦!  + Vp-  64 

! 

(45)  } 

"t 

^ ♦*w-«3 

(46)  J 

. ! 

then 

> 

! 

' i 
\ 

■ i 

• l 

Ff.  • 0 

Cf  > 0 

1 

\ 

1 

j 

1 

■ i 

■ -K.topGfK3 

0<G ,<  h. 

(47)  ; 

* 

l 

i 

* ‘“.top  «W  K3  + K.t.p  K3  h. 

Qf>h, 

1 i 

j 

Fb„-  0 

G„>  0 

- i 

! 

■ -K.topabK3 

0«Ob<  h. 

(48) 

. * 

'.1 

* ’".top  «Vh.>  K3  + K.top  K3  h. 

°b>h. 

■ i 

Where  K#t  is  the  stop  stiffness  factor  for  bogie- stop 
and  hgtQp  is  the  stop  stiffness  factor  for  bogie-hull 
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APPENDIX  C 


VEHICLE  RESPONSE  PLOTS 


A complete  set  of  PSDs  of  vehicle  response  obtained  from  the 
M-809  truck  simulation  are  ehown  in  Figures  C-l  to  C-10.  The 
following  symbols  are  keyed  to  the  various  curves. 

Squaro  Point  Contact  Tire  Model 

Circle  Rigid  Tread  Band  Tire  Model 

Triangle  Fixed  Footprint  Tire  Model 

Cross  Adaptive  Footprint  Tire  Model 


Power  Spectral  Densities  of  Vehicle  Force.  Motion  and 
Acceleration  for  M-809  Truck  (Speed  - 18  mph.  Terrain  Roughness  - lnrmi) 


Cg  Acceleration 

C-l 

Front  Axle  Acceleration 

C-2 

Rear  Axle  Acceleration 

C-3 

Front  Arle  Displacement 

C-4 

Middle  Axle  Displacement 

C -5 

Rear  Axle  Displacement 

C-6 

Vertical  Tire  Force  (Front) 

C-7 

Vertical  Tire  Force  (Rear) 

C-8 

Fore  -and -aft  Tire  Force  (Front) 

C-9 

Fore-and-aft  Tire  Force  (Rear) 

C-10 

PSD  (FT  /SEC  /RAD)  : 


CG  flCCLN 


□ Point  Contact  Tire  Model 
O Rigid  Tread  Band  Tire  Model 
A Fixed  Footprint  Tire  Model 
•f  Adaptive  Footprint  Tire  Model 


Frequency-,  hs 
10 


REDUCED  FREOUENCY  - RflD/FT 


Figure  C-l.  Simulated  Response  - CG  Acceleration 


9 


/SEC  /RAD) 


/SEC  /RAD) 


REAR  AXLE  flCCLN 


□ Point  Contact  Tire  Model 
O Rigid  Tread  Band  Tire  Model 
A Fixed  Footprint  Tire  Modol 
4.  Adaptive  Footprint  Tire  Model 


Frequency,  he 


REDUCED  FREQUENCY  - RfiD/FT 


Figure  C-3,  Simulated  Retponae  - Rear  Wheel  Acceleration 


PSD  (FT  /RAD) 


FRONT  RXLC  DISPLN 


Figure  C-4.  Simulated  Response  - Front  Wheel  Displacement 
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MIDDLE  AXLE  DISPLN 


□ Point  Contact  Tire  Model 
O Rigid  Tread  Band  Tire  Model 
A Fixed  Footprint  Tire  Model 
4.  Adaptive  Footprint  Tire  Model 


Frequency,  ha 


PSD  (FT/RAD) 


REAR  AXLE  D1SPLN 


□ Point  Contact  Tire  Model 
O Rigid  Tread  Band  Tire  Model 
A Fixed  Footprint  Tire  Model 
4.  Adaptive  Footprint  Tire  Model 


Frequency,  he 


REDUCED  FREQUENCY  - RRD/FT 


Figure  C-6.  Simulated  Reeponee  - Rear  Wheel  Displacement 


PSD  (LB  FT  /RAD) 


FRONT  TIRE  FORCE  (vertical) 


□ Point  Contact  Tire  Model 
O Rigid  Tread  Band  Tire  Model 
& Fixed  Footprint  Tire  Model 
4.  Adaptive  Footprint  Tire  Model 


Frequency,  hz 


REDUCED  FREQUENCY  - RRQ/FT 


Figure  C-7.  Simulated  Response  - Vertical  Tire  Force  (Front) 


PSD  {LB  FT/RAD) 


REAR  TIRE  FORCE  (vertical) 


□ Point  Contact  Tire  Model 
O Rigid  Thread  Band  Tire  Model 
A Fixed  Footprint  Tire  Model 
4 Adaptive  Footprint  Tire  Model 


Frequency,  hz 


X 1“ 

10 


— i — i -r~rm  | 1 1 — rrnrn 1 r — » Trrrri 

10°  10'  n 

REDUCED  FREQUENCY  ~ EDO /FT 


Figure  C-8,  Simulated  Response  - Vertical  Tire  Force  (Rear) 


RERR  TIRE  FORCE  (FORE-AND-AFT) 


□ Point  Contact  Tire  Model 
O Rigid  Tread  Band  Tire  Model 
A Fixed  Footprint  Tire  Model 
+ Adaptive  Footprint  Tire  Model 
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Figure  C-10.  Simulated  Responae  • Fore-and-Aft  Tire  Force  (Rear) 
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Rtcorded  simulation  control  scheme  deecribed  in  this  report  is  s 
method  of  providing  realistic  inputs  to  a vehicle  shaker.  The  scheme 
involves  generating  axle  displacement  records  by  simulating  motion  of 
the  vehicle  operation  on  the  specified  terrain,  storing  the  records  in 
memory  of  a control  ayatem  and  using  them  to  provide  input  signals  to 
the  vehicle  theker.  , , (Continued) 
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